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ABSTRACT 


Analytical methods have been developed to predict the wave motion in 
baffled combustion chambers. Principally, the method involves solution 
of the wave equation through use of an iterative/variational method. 
Calculated results are in good agreement with available data from bench 
scale acoustic models. Convergence of series expressions which occur 
in the analysis and of the iteration scheme was demonstrated numerically. 
Calculated results with gain/loss end- wall boundary conditions do not 
agree with observed engine stability, which indicates that other gain/ 
loss representations must be used. The principal contribution from 
this study to aid the design of baffles is the ability to accurately 
predict instability frequencies and oscillatory pressure and velocity 
distributions for baffled chambers. 
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NOMENCLATURE* 

= sound velocity, in. /sec 
= Green's function 
= to/c, inch ^ 

= (k 2 - q 2 ir 2 /w 2 )^ 2 , inch 1 

= (k 2 - m 2 ^ 2 /^ 2 ) 1 / 2 , inch" 1 

= baffle length, inches 

= main chamber length, inches 

= (£ + L), total chamber length, inches 

= main chamber summation index, positive integral values 

= steady flow Mach number 

= baffle compartment summation index, positive integral value 
= unit normal vector directed outward 
= pressure, lb f /in? 

= summation index, positive integral values 
= position vector, inches 
= radial coordinate, inches 
= chamber radius, inches 
= surface area, sq in. 

= width of baffle compartment, inches 
= width of main chamber, inches 
= longitudinal position coordinate, inches 


*A separate table of nomenclature is included in appendices C and D. 
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a) 


= transverse position coordinate , inches 

= specific acoustic admittance of injector end (dimensionless) 

= specific acoustic admittance of nozzle end (dimensionless) 

= axial coordinate , inches 

= damping coefficient , imaginary part of complex angular frequency, 
seconds"*' 1 ' 

= 1 if v = o, 2 if v / 0 

= baffle compartment index, positive integral values 
= normal pressure gradient, lb^/in, 3 

2 

= time averaged gas density, lb /in, 

= coW/c or cjR/c 

= angular frequency, radians/sec 
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to main chamber 

to baffle compartments 

to radial summation 

to interface plane 
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to source coordinates for Green ! s functions 
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SUMMARY 


Analytical methods have been developed to predict the wave motion in 
baffled combustion chambers. These methods may be usefully applied to aid 
the design of injector face baffles. 

An analytical method has been developed by employing variational and 
iterational techniques to solve the wave equation for baffled chambers. 

For this purpose the wave equation and boundary conditions were converted 
to an integral equation which was solved by these approximate means. The 
analytical method has been applied to both two-dimensional and cylindrical 
chambers with no steady flow and with relatively simple baffle configura- 
tions. Results from this analysis provide a good estimate of the instability 
frequencies and oscillatory pressure and velocity distributions in similar 
combustion chambers. Application of the analysis to more complicated 
baffle configurations appears straight-forward. The analysis was partially 
extended to include uniform steady flow in the chamber, but no numerical 
results were obtained. Nonetheless, no substantial difficulties associated 
with this extension were found. 

Analytical results for closed chambers were found to be in excellent agree- 
ment with available data from bench-scale acoustic models. In addition, 
convergence of the several series-type expressions which appear in the 
analysis and, also, of the iteration scheme was investigated, retailed 
calculations were used to demonstrate satisfactory convergence of the 
several series and of the iteration scheme. 
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Stability calculations were made for a number of chamber/baffle configurations 
with gain/loss end-wall boundary conditions. Admittance -type boundary con- 
ditions were used with signs chosen to represent a loss at the nozzle end and 
a gain at the injection end of the chamber. Results from each of the cases 
analyzed showed that, as long as the injector had a gain-type admittance, 
the predicted stability of the chamber was degraded by the presence 
of baffles. Because this result is obviously incompatible with observed 
engine stability, this simple representation of gains and losses must be in- 
adequate. This failure to properly predict observed stability trends may be 
due to effects of velocity-coupled combustion driving or of wake-type drag 
losses on the baffles, neither of which are included in the current analytical 
formulation. 

The principal contribution from this study to aid the design of baffles is 
the ability to accurately predict instability frequencies and oscillatory 
pressure and velocity distributions in baffled chambers. 
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INTRODUCTION 


The intent of this report is to describe recently developed analytical 
methods and results which may be used to aid the design and development of 
injector face baffles. These results and methods have been obtained from 
an analytical study to investigate the stabilizing effects of baffles on 
acoustic modes of combustion instability. Many of the analytical details 
and supporting studies obtained from this investigation will be described 
only briefly so that emphasis may be placed on information useful for baffle 
design and development. 

Injector face baffles are often required to prevent acoustic modes of combus- 
tion instability in rocket engines. These baffles are currently designed on 
a largely empirical basis. Some attempts have been made to systematize avail- 
able experimental results (for example, see Ref. 1); other attempts have been 
made to use the Priem-type instability model to guide baffle design (Ref. 2). 
These approaches still provide only rough design guidelines, Analytical 
models are now available which will predict the stabilizing effects of baffles 
and, thereby, provide a theoretical basis for baffle design. 

The acoustic modes of instability for unbaffled engines are currently approxi- 
mated by the corresponding modes of a similarly shaped closed acoustic cavity, 
e.g., a cylinder. Employing methods described herein, the acoustic modes of 
closed but baffled chambers can also be calculated and used for similar purposes. 

Calculation of the acoustic modes of closed baffled chambers was only an interim 
goal of the program on which this report is based. Nonetheless, this capability 
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can be used to aid baffle development by providing approximate frequencies, 
and oscillatory pressure and velocity distributions for the instability modes 
that can occur in the baffled engine. Thus, the general acoustic mode charac- 
teristics can be understood even though stability of the engine cannot be 
predicted. This analysis shows that baffled chambers may exhibit well-defined 
acoustic modes similar to the corresponding modes of unbaffled chambers. This 
conclusion is supported by available data from bench-scale acoustic models and 
hot firings. Thus, the normal acoustic modes are not el imin ated by the intro- 
duction of baffles, they are only disguised. 

ANALYTICAL APPROACH 

The analytical approach being used has been generally described in Ref. 3. 
Without steady flow the analysis concerns approximate solution of the wave 
equation, which simply represents a composite of the linearized fluid dynamic 
equations obtained through the assumption of small variations from the mean 
(time average) value. With uniform steady flow, additional terms arise so 
that an inhomogeneous wave equation is obtained but, nevertheless, the general 
approach remains unchanged. 

Solution of the wave equation for a baffled chamber is complicated by the fact 
that the standard separation-of -variables technique cannot be used because of 
the boundary shape. For similar problems, Morse (Ref. 4, pg. 1039) suggests 
the use of an integral formulation which may be subsequently solved by approxi- 
mate means. This suggestion has been followed in this study. 

The wave equation and boundary conditions were rewritten as an integral equation 
which, in turn., was solved by a combined variational and iterational method. 
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This approach was found to work very well. The approach has been applied to 
both two-dimensional and cylindrical baffled chambers and has been extended 
to include uniform steady flow in the two-dimensional case. Further, the 
effects of combustion driving and nozzle losses have been simulated to 
some extent through incorporation of gain/loss-type boundary conditions. 

Integral Formulation 

The wave equation and accompanying boundary conditions were converted to an. 
integral equation through the use of Green's functions, as described by 
Morse (Ref. 5, pg. 321). The same conversion can be accomplished for this 
case without the introduction of Green's functions if a properly formulated 
sum of separated- type solutions is used (described in Appendix A). None- 
theless, the use of Green's functions greatly facilitates the analysis 
and allows it to be done more systematically. 

The Helmholtz equation (which is the wave equation for a harmonic time 
dependence), i.e., 

V 2 p + k 2 p =0 (l) 

may be rewritten as 

P (r) = f G(r|r ) N • V Q p(r )dS (2) 

Jq ° u o o 

where G(r|r ) is a Green's function, which satisfies either the same boundary 
conditions as the pressure (p), or a zero-gradient boundary condition. In 
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addition. Green's function satisfies the differential equation 

V 2 G + k 2 G = -6(r-r ) (3) 

where 5i r-r ) is a Dirac delta function. Expressions for the Green's func- 

o 

tion may he obtained in several ways, such as described by Morse (Ref. 4, 

pp 791 to 834). 

The integral expression for pressure is used with separate Green's functions 
written for each baffle compartment and also for the main chamber. Each of 
these Green's functions is zero outside of the compartment to which it 
applies. However, the oscillatory pressure and normal component of velocity 
must be continuous across the conceptual interface between each region. 
Therefore, at this interface, 

P (r_) = f G (r |r )£(r )dS = - f G, (rlr )i(r_)dS = p, (r ) (4) 

a s' J a ' s o' 0 ' o J bp' s 1 o o' o ^bp v s / ' 

S 0 

where G„(r|r„) is the Green's function for the main chamber, G, (rlr ) 

S- 'O DfJ. . 0 

is the Green's function for the p^ 1 baffle compartment, and, to simplify the 

notation, 

§(?o) = N • p a (^) = - H*V 0 p b (^) (5) 

In addition, the normal pressure gradient, |, must satisfy the integral 

equation 

Ut 5 ) =/S-v 0 G a (? s |^)5(^)ds 

( 6 ) 
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Equation 6 was obtained directly from Eq. 2 by differentiation. 

Simultaneous solution of Eq. 4 and 6, or their equivalents, to give the 
allowed frequencies of oscillation and the normal pressure gradient is the 
pivotal portion of the analysis. Both the Green's functions and i depend 
on frequency. With this information, the oscillatory pressure and velocity 
at any point in the chamber may be calculated from the integral expression 
for pressure, Eq. 2, by integration. The velocity components may be obtained 
from the gradient of the pressure. A variational-iterational technique has 
been developed to solve Eq. 4 and 6. 

This approximate solution technique results from a direct combination of a 
variational technique, which allows the "best" form of an approximate solu- 
tion to be selected, and an iteration technique, which allows an initially 
selected approximate solution to be systematically improved (also see 
Ref. 3). The variational procedure results in replacement of Eq. 4 by a 
characteristic equation; the procedure used is described by Morse and 
Ingard (Ref. 5, pg. 68o). Employing the variational function developed 
by them, with a slight generalization for multiple compartments, a charac- 
teristic equation was obtained. 



Where an approximate normal velocity of the form u = A£, has been used, the 
value of the amplitude (A) being optimized by the method. Thus, at this 
level of approximation, by assuming a reasonable estimate for the normal 
gradient, the frequency may be calculated from Eq. J. Employing the 
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Integral expression for pressure (Eq. 2) Eq. 7 may be rewritten as 

£ (P a -P b )<*S * o (8) 


L 


If the exact expressions for £ is obtained, then Eq. 8 is satisfied identically 

because p = p . However, as it has been used, the variational procedure 
a b 

does not indicate how to estimate £; rather, it indicates that the best 
estimate of the allowed frequencies corresponding to a particular estimate 
of £ is obtained by satisfying Eq. 8. The iteration procedure is used to 
obtain an arbitrarily good estimate of £ . Nonetheless, continuity of 
pressure at the interface is satisfied only in the average sense defined by 


m P, 

JuhJw w 1 * 


Tine iteration procedure has been set up so that an initial estimate 


>f £ may be improved with iterations such that in the limit p = p . 

a b 


The iteration procedure consists of assuming a pressure distribution for 
p , calculating the corresponding £ and then p from that £ (by employing 


a 


the integral expression for pressure). A new estimate is obtained for p^ 

by equating it to the newly calculated p . All of this is possible because 

a 

the expressions for p , p and £ are expressed as series of orthogonal 

a t, 

functions. 


This procedure leads to series-type algebraic expressions for the character- 
istic equation and the pressures. These axe easily solved by numerical 
means. By such means, convergence of the iteration scheme has been demon- 
strated, i.e., p does indeed approach p, with iterations. However, other 
variations of this scheme did not converge, as will be discussed subsequently. 


8 



Gain/Loss- Type Boundary Conditions 


Nonzero , admittance -type boundary conditions can be added at each end of 
the chamber without difficulty. By defining Green’s functions, which 
satisfy the same boundary conditions at the closed ends of the chamber as 
those satisfied by the pressure, the foregoing equations may be used without 
change . 


Uniform Steady Flow 


If a uniform steady flow is assumed to exist in the chamber, the analysis 
must be modified to some extend. The linearized fluid dynamic equations 
for small amplitude oscillations may be combined to give an inhomogeneous 
wave equation. 


V 2 p + k 2 p - 2jkM 9£. 

° dz 


2 

d P 
o o 
dz * 


M 2 


(9) 


This equation with the boundary conditions can be converted to an Integral 
equation as before and a parallel analysis may be applied. During the 
course of the program, the characteristic equation and iteration equations 
were developed but no numerical results have been obtained. 
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' ANALYTICAL RESULTS 


These analytical methods have "been used to investigate the oscillatory 
characteristics of both two-dimensional and cylindrical baffled chambers 
with either rigid end- walls or gain/loss end- wall boundary conditions. 
Further^ the convergence properties of the iteration scheme and of the 
several series expressions which arise in the analysis have been care- 
fully examined. Finally, the analysis has been extended to include the 
influence of uniform steady flow but no numerical results have been obtained. 

The utility of these methods are dependent upon the convergence properties 
of the expressions. To be useful, acceptable convergence must occur with 
reasonable computing requirements since the equations are solved numerically. 
Generally, the convergence properties were found to be good and the com- 
puting requirements (computing time or storage) were not considered excessive 
for any case. These factors will be discussed more completely later in this 
report. 

Rigid Walled Chambers 

The oscillatory characteristics of several baffled chambers with rigid walls 
(no gains or losses) and no steady flow were examined in a series of calcula- 
tions, A number of calculations were made for two dimensional chambers con- 
taining one or more zero-width baffles of equal length and spacing. A number 
of calculations were also made for cylindrical chambers containing three 
radially directed zero-width baffles of equal length and spacing. Results 
from the latter calculations are probably of the greatest practical interest. 
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These calculations were done for the chamber configurations shown in Fig. 1 
and 2. These also serve to indicate the dimensional notation used in the 
analyses. The relatively complicated expressions used to make these calcu- 
lations are given in Appendices C and D along with listings of the corre- 
sponding computer programs. The calculations were restricted to the first 
transverse and the first tangential modes. 

Typical pressure and velocity profiles for the two-dimensional case are shown 
in Fig. 3* The calculations clearly indicate a discontinuity (singularity) 
in the axial component of velocity at the baffle tip. In the physical situa- 
tion, this discontinuity would be replaced by a steep gradient region. 

Clearly a strong flow around the baffle tip is indicated. Also shown is the 
degree to which a pressure match has been achieved with the indicated ten 
iterations. Note that the pressure distributions are roughly similar to the 
sinusoidal distribution obtained without a baffle. The pressure amplitude 
was greatest at the injector plane and diminished with length, as shown in 
Fig. 4. 

Typical pressure and velocity distribution for the cylindrical case are shown 
in Fig. 5 and 6. Again, a discontinuity in the axial component of velocity at 
Ihe baffle tip is indicated. The oscillatory character of this velocity dis- 
tribution on either side of the baffle would probably be diminished if a 
greater number of terms in the series were retained, i.e. a more nearly 
monotonic distribution is expected. Note that the pressure match is reason- 
ably good. 
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Figure 1. Configuration and Notation Used in 
Two-Dimensional Chamber Analysis 
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Axial Velocity - Baffle Compartment 
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Transverse Profiles at the Interface 
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Interface Pressure and Velocity Profiles in Radial Direction (0 - 0+) 





Cylindrical Chamber With Three Equally 
Spaced Radial Baffles 
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The calculated frequencies were found to be depressed from the unbaffled values 
as the baffle -length/ chamber -length ratio was increased or as the chamber- 
width (radius) /chamber -length ratio was increased. These effects are illustrated 
by the frequency results shown in Tables land 2 and in Fig. 7* 

It should be noted that solution of the characteristic equation to obtain the 
frequencies (the roots or zeros of this equation) yields more than one frequency 
in the range of interest. This effect is illustrated in Fig. 8 which shows 
the variation of the characteristic equation with frequency. Two roots or 
allowed frequencies are shown, only one of which corresponds to the first trans- 
verse mode. The second root, usually the lower one, is discarded because it 
corresponds to a poor pressure match at the interface. This discarded root 
say represent a poor estimate of a physically realizable frequency or may 
simply be an extraneous root. 

The multiple roots arise simply because there are different combinations of 
series coefficients (these are frequency dependent) which will satisfy the 
characteristic equation. Generally, this effect has not been troublesome be- 
cause once the appropriate root is located its variation can be followed when 
various parameters are changed. The criteria for selecting the roots to be 
utilized should always be the degree of pressure match. 

By and large, all of the foregoing discussion has concerned calculations with 
the normal gradient, i , represented as a series of compartment eigenfunctions. 
Other calculations have been made with | expressed as a series of the main 
chamber eigenfunctions and, also, with alternating .expressions. None of the 
latter calculations produced results as good as those described, a poor pressure 
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TABLE 1. CALCULATED EIGENVALUES FOR TWO-DIMENSIONAL CHAMBERS 
WITH ZERO ADMITTANCE BOUNDARY CONDITIONS* 


Objectives 

V w 

c/i T 

00 W 

C: 

Effect of ■fi/L 

1.5 

0.1 

3.05217 

T 

1.5 

0.2667 

2.3742 


1.5 

0.4333 

1.76S41 


1.5 

0 . 5667 

1.44435 

Effect of L^/w 

1.0 

0.25 

2.783564 : 

1.5 

0.25 

2.451465 


2.0 

0.25 

2.103642 | 

Effect of No. of Baffles 




1 - Baffle 

1.5 

0.2667 

2.3742 

4 - Baffles 

1.5 

0.2667 

2.24889 


*q = 10; m = 15; 10 iterations 


TABLE 2. CALCULATED EIGENVALUES FOR CYLINDRICAL CHAMBERS 
WITH ZERO ADMITTANCE BOUNDARY CONDITIONS* 


Objectives 

V R 

C/Ivp 

coW 

c ! 

Effect of C/l 

2.8 

0.00 

1.84118 | 

' T 

2.8 

0.10 

1.80062 ! 


2.8 

0.25 

1.44294 1 


2.8 

0.50 

O.9038I j 

Effect of L^/R 

1.0 

0.25 

1.76668 j 

2.8 

0.25 

1.44294 | 


5.0 

0.25 

0.96102 


*m 0 = 5, n 0 = 5, m r = 3, = 3; 10 iterations 
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Figure 7. Predicted Frequency Dependence on Baffle Length 
From the Iterative Approximation 
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match being obtained. This difference in results stems frcm the velocity 
singularity that occurs at the baffle tip, as shown in Fig. 3. The compart- 
ment series can readily describe this singularity because a separate expres- 
sion is used for each compartment and it need not be continuous between 
compartments. The main chamber series, on the other hand, must be continuous 
at the baffle tip. Consequently, many more terms are required in the main 
chamber series than in the compartment series to approximate the velocity 
singularity. It is entirely possible that, if a sufficient number of terms 
is retained in the main chamber formulation, comparable pressure matches could 
be obtained, but no clear indication of this has been observed. Thus, it may 
also fail. 

Convergence 

A number of calculations were made to investigate the convergence properties of 
the analytical method . Generally this investigation was directed toward the 
two-dimensional rigid walled chamber. The results are believed generally 
applicable to the cylindrical case as well. 

Two kinds of convergence considerations are required, because the set of 
equations developed to describe the baffled chamber contains several infinite 
series as well as an iterated approximation. The infinite series arise from 
the series representations for the oscillatory pressure and velocity in the 
main chamber and separate but similar representations for the baffle compartments. 
These series arise irrespective of the iterated approximation. Therefore, 
independent of whether or not the iteration scheme converges, a sufficient 
number of terms must he retained in each of the infinite series to adequately 
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approximate the entire series. Thus, for a two-dimensional chamber convergence 
is required of the series in m and q (series indices for the main chamber and 
compartment expansions, respectively). Once the required numbers of terms 
(m and q) were determined, the question of iterative convergence could be 
investigated. 

Generally it was found that relatively small numbers of iterations and terms 
(i.e. terms retained in the various series expressions) were required before 
the calculated oscillatory frequencies ceased to change as more terms and 
iterations were added. Somewhat larger numbers of terms and iterations were 
required before the calculated pressure distributions ceased to change. The 
calculated velocity distributions, especially near the baffle tips, required 
the largest numbers of terms and iterations. 

Initially convergence was investigated by examining the variation of the 
pressure profile on the compartment side of the chamber/ compartment interface 
as the number of terms in q was increased but with fixed numbers of terms in 
m and iterations. When this profile became insensitive to further increase, 
convergence of the series in q was considered to be achieved. A similar 
procedure was then used to determine the required number of terras in m, 
this time with the number of terms in q fixed at the level required for con- 
vergence of that series. Finally, with an adequate number of terms in m and 
q, the number of iterations was increased. 

For a single baffle in a two-dimensional chamber with a baffle-length- to- 
chamber-length ratio of 0.1 and a chamber length- to- width ratio of 1.5, 
results from the indicated set of convergence calculations is shown in Fig. 9a 
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Oscillatory Pressure Compartment 









through 9d. In Fig. 9a compartment- side pressure profiles are shown for 
various q, i.e., number of terms in q. Corresponding calculations are shown 
in Fig. where the number of terms in m was varied. Clearly convergence 
is occurring and approximately 10 terms in q and 15 terms in m are adequate 
for this case. It should be noted that for the first transverse mode all 
of the even order terms in m are zero, which is required for proper symmetry 
of the pressure profile. 

Figures 9c and 9d show comparative pressure profiles on each side of the 
chamber/compartment interface for various numbers of iterations. Again, it 
is clear that convergence is occurring. Error calculations show that the 
error decreases in a monotonic fashion. Errors are even smaller with more 
than one baffle. 

Somewhat larger numbers of terms are required to describe the velocity pro- 
file near the baffle tip. In Fig. 10 are shown calculated profiles from 
the truncated series with q = 10 and m = 15 . Clearly, the series for the 
main chamber velocity is not adequately approximated with 15 terms. How- 
ever, Fig. 3 demonstrates that an excellent velocity match is achieved if 
m is increased to m = 45. Figures 11 and 12 show calculated results for 
10 and 20 iterations, respectively, and with q = 23 and m = 45 . The velocity 
agreement would improve if the ratio of m to q was increased. 

Adequate convergence was considered to have been demonstrated by these cal- 
culations. Additional convergence considerations are described in Appendix B. 
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1 Baffle 

Figure 11, Calculated Pressure and Axial Velocity At interface in Two-Dimensional 




(sisog iCjBaq.fqay) Xq-pooiaA ao ajnssaOjj 


28 


Two-Dimensional Chamber 




Gain/ Loss Boundary Conditions 


The effects of pressure-coupled, combustion driving and nozzle losses were 
partially simulated by employing nonzero acoustic admittance values at each 
end of the chamber. A gain-type condition was imposed at the injector end 
and a loss-type condition at the nozzle end. 

The boundary condition was written in terms of a specific acoustic admittance: 

y = pc jj. * u 
P 

with y T being the injector admittance and y being the nozzle admittance. 
x H 

With such boundary conditions, the characteristic equation becomes complex 
and solution of it yields complex frequencies or eigenvalues. The real part 
of that eigenvalue corresponds to the angular frequency of oscillation; the 
imaginary part corresponds to an exponential growth or decay coefficient. 

Thus, whether stability or instability is predicted depends on the sign of 
this imaginary part of the eigenvalue. 

A series of calculations were made for both two-dimensional and cylindrical 
chambers. The results are shown in Tables 3 and 4. In every case, the 
baffles were found to degrade stability (as long as the real part of yj was 
negative) rather than improve it as found in practice. This implies that the 
stability model is inadequate, possibly because velocity (transverse) coupling 
with the combustion is more important than the pressure sensitive case that 
was simulated. Another possible explanation is that the baffles contribute 
nonlinear losses due to wake effects which were not included in this analysis,. 


29 



TABLE 3, DAMPING CALCULATIONS FOR THE NONZERO ADMITTANCE CASES* 


Objectives 

Number of 
.Baffles 

w/l t 

! /h 

y N 

y i 

co W 
c 

a W ; 

c 

Effect of 

1 

1.0 

0.25 

0.06 

-0.06 

2.722 

-0.0773 

number of 

2 

1.0 

0.25 

0.06 

-0.06 

2.706 

-0.0527 

baffles 

3 

1.0 

0.25 

0.0 6 

-0.06 

2.654 

-0.0557 

' Effect of 

3 

1.0 

0.10 

0.06 

-0.0 6 

3.059 

-0.0153 

1 baffle 

3 

1.0 

0.25 

0.0 6 

-0.06 

2.654 

-0.0557 

1 length 

3 

1.0 

0.50 

0.06 

-0.06 

1.925 

-0.0599 

Effect of 

3 

1.0 

0.25 

0.06 

-0.20 

2.680 

-0.2170 

injector 

3 

1.0 

0.25 

0.06 

-0.0 6 

2.654 

-0.0557 

admittance 

3 

1.0 

0.25 

0.06 

-0.02 

2.652 

-0.0106 

Effect of 

3 

1.0 

0.25 

0.2 

-0.06 

2.665 

-0.0305 

nozzle 

3 

1.0 

0.25 

0.06 

-0.06 

2.654 

-0.0557 

admittance 

3 

1.0 

0.25 

0.02 

-0.06 

2.653 

-0.0634 

Effect of 

3 

2.0 

0.25 

0.06 

-0.06 

2.896 

-O.OI76 

chamber 

3 

1.0 

0.25 

0.06 

-0.06 

2.654 

-0.0557 

width. -to« 
iengtfa 

3 

0.50 

0.25 

0.06 

-0.06 

1.996 

-0.0702 


*A11 calculations were made with 15 terms in m, 10 terms in q.j, and 

9 iterations. 


TABLE 4 . CALCULATED EIGENVALUES FOR CYLINDRICAL CHAMBERS 
WITH FINITE ADMITTANCE BOUNDARY CONDITIONS* 


Objectives 

V R 


y N 

y i 

Co r 
w 

c 

ar 

w 

c 

Effect of fi/L 

2.8 

0.10 

0.06 

-0.06 

1.80472 

-O.OI665 

1 

2.8 

0.25 

0.06 

-0.06 

1.H383 

-0.04459 

1 

2.8 

0.50 

0.06 

-0.06 

0.90390 

-0.03392 

Effect of ly'R 

1.0 

0.25 

0.06 

-0.06 

’ 1.77020 

-0.00729 

2.8 

0.25 

0.06 

-0.06 

1.44383 

-0.04459 


5.0 

0.25 

0.06 

-0.06 

0.98895 

-0.03405 

Effect of 

2.8 

0.25 

0.00 

0.00 

1.44294 

«•»«. i 

1 End Admittance 

2.8 

0.25 

0.06 

-0.06 

1.44383 

-0.04459 

[ 

2.8 

0.25 

0.20 

-0.20 

1.45344 

-0.15259 


*tn„ = 5 > m = 3> n =* 3> 10 iterations 

6 O’ r r 
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The results from stability limit calculations made for the two-dimensional 
case with gain/loss boundary conditions were similar to those 
reported previously in Ref. 3j no region was found in which a beneficial 
effect of baffles was indicated. 

Computing Time Requirements 

The utility of the foregoing analytical methods depends in part on the com- 
puting requirements in terms of computing time and storage. Therefore, it 
is of interest to determine these requirements for a "typical" case. Un- 
fortunately, a typical case is very difficult to specify because of the 
variability in the numbers of terms and iterations that may be dieemed 
necessary. Moreover, occasionally considerable difficulty is encountered 
in finding the roots of the characteristic equation. Projecting the require- 
ments of another computer system is even more difficult. Nonetheless, an 
attempt has been made to estimate the computing time requirements. 

All calculations have been made on a G.E. Model 440 time sharing computer 
system. For "typical" cases the rough maximum requirements are shown below. 

Two-dimensional 

Rigid wall case < 150 sec or 

< $0.42 

Nonzero admittance case < 900 sec or 

< $2. 50 


Cylindrical 

< 450 sec or 

< $ 1.25 

< 1500 sec or 

< $ 4. 20 
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'Hie dollar figures shown are based on a nominal rate of $10 per hour of 
computing time. The salary of the engineer operating the computer is not 

included® 

If these same cases were to non on an IBM-360 Model 65 computer in a batch 
processing rather than time sharing mode, the computing times would be l/2 
(very roughly) of the times shown above. This estimate includes an allowance 
for requiring more trial cases because of the inability to interact. directly 
with the computer. The computer logic being used has, of course, been 
selected for a time sharing mode. Therefore, some improvements are possible 
for use with a batch processing computer. 

Clearly the computing costs are well within an acceptable range. 


CONCLUSIONS 

A satisfactory analytical method (variational-iterational technique) has been 
developed to analyze the wave motion in baffled combustion chambers. This 
method may be used to accurately predict the instability frequencies a n d , the 
oscillatory pressure and velocity distributions in baffled chambers. A 
number of relatively simple configurations have been analyzed thus far and 
application of the method to other baffle configurations appears straight- 
forward, although the details of the analysis are complicated. 

Stability calculations have been made for a number of chamber/baffle configura- 
tions with gain/loss end-wall boundary conditions, which were used to simulate 
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combustion gains and chamber losses. In each case analyzed (with gain type 
admittance at the injector end), the predicted stability of the chamber was 
degraded by the introduction of baffles. Because this result is in direct 
conflict with observed engine stability, this simple analytical representa- 
tion of chamber gains and losses must be inadequate. This failure to 
properly predict observed stability trends may be due to velocity-coupled 
combustion driving effects or to wake -type drag losses on the baffles, 
neither of which are included in the current analytical formulation. 

Nevertheless, the ability to predict instability frequencies, and oscillatory 
velocity and pressure distributions in baffled chambers may be used to 
advantage for the design of baffles. 
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APPENDIX A 


ALTERNATE DERIVATION OF INTEGRAL EQUATION 


The integral equation for pressure can he derived from the wave equation by 
employing a sum of separated solutions rather than starting with a Green's 
function. Essentially this approach is an alternate derivation or definition 
of the surface Green's function. The method is described for the two-dimen- 
sional case but the cylindrical closely parallels it. 


A two-dimensional chamber is considered with rigid walls on three sides 
but an arbitrary boundary condition at x = L. The separation of variables 
technique is used without applying the boundary condition at x = L. The 
separated solutions may be written as 


f n (*,y) ■ XX. eos kx (A 

2 . 2-2 2 l /2 

where = (k -mfr/w)'. Note that only one index, m, appears. Because 
the individual separated solutions cannot satisfy the arbitrary boundary 
condition at x = L, a sum of these solutions is tried, i.e.. 


/ 


p = 



mrry 

W 


cos k x 
m 


(A-2) 


where a^ is an unspecified coefficient. The pressure gradient at x = L may be 
obtained by differentiating Eq. A-2, 


9p 

8x 


x = 


L 



^ . . . T IU/J.Y 

a k sm k L cos —rr~ 

m m m W 


(A-3) 


A-l 



Th 


efficients may be obtained in terms of this gradient from orthogonality, 

r W Or 


3x 
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cos 


X = L 


m7r y 0 

~w 


dy 
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W 

— k sin k L 
e m m 
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2 m7ry 
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m = 0 
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Thus, the pressure at any point in the chamber may be expressed in terms of 


the gradient at x = L, i.e., Eq. A-2 becomes 

JW 


, m7ry 

cos k x cos ~- 
m W 


"m 


/. 


9P 

ax 


x = L 


COS mV y° <}y 
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m 


sin k^L 


(A-5) 


By noting that the integration in Eq. A-5 is over the variable y^ and that a 
series may be integrated term by term, Eq. A-5 may be written as 


3x 


T G dy D 

x = L O 


(a-6) 


A-2 



where, by definition, 


G = 


I 


' m 


W 


cos 


cos 

m 


m7ry 


cos 


mny r 


W 


k sin k L 
ra m 


(A-7) 


Equations A-6 and A-7 are identical with those obtained by starting with a 
Green's function (see Ref. 3 on Appendix C). 



APPENDIX B 


ADDITIONAL CONVERGENCE CONSIDERATIONS 


It is important to define the degree to which a continuous oscillatory velo- 
city distribution is obtained with the analytical method being used. This 
question relates closely to convergence of the series representations for 
pressures and velocities, but not (directly, at least) to convergence of the 
iteration scheme. In general, each of these series should be regarded as 
being of infinite extent but are approximated by truncated series for 
numerical purposes. 

The analysis is such that if the normal gradient distribution at the chamber/ 
compartment interface ( % ) and the frequency (or k) are given, then the pres- 
sures and pressure gradients (or velocities) can be subsequently calculated 
for any point in the chamber or compartment from the integral expressions. 

For example, the axial gradient of pressure in the main chamber may be written 
(for the two-dimensional case) as 


d p a V Jtn 
9x Z W 


sm 


k x 
m 


SIB 


k L 
m 


cos 




(B-l) 


If x = L, this expression is simply a Fourier series expansion for % . A 
similar result is obtained for the baffle compartments, i.e., 

dy 

o 

(B~2 


2c 

dx 
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£ sin k (L, - x) 

h g t / 
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sin k S. 
< 3 . 


cos 


q rr (y - ;uw) 
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J ' % cos 


q*r(y -mw) 
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Again, for x = L, the expression becomes a Fourier series expression for £ . 


Because for x= L each of these expressions is a Fourier series for the same 
thing, then the normal component of velocity is continuous, at least if the 
series converge. Any deviation from continuous character must arise from 
truncation of the series for the calculations. This conclusion is valid for 
any level of iteration or approximation for £ (and the solution); however, 
the pressure is not continuous unless the exact solution is obtained. 

The variational-iterational technique is the means being used to calculate 
approximations for £ and k. With the compartment approximation, the Fourier 
series for I in terms of the compartment-side eigenfunctions is used, which 
may be denoted as £ . The eigenvalue (frequency), in turn, is calculated 
from the characteristic equation, 

jt't - V h )^-0 (B 

The iteration is used to generate approximate expressions for p (l), 
and p^(l); the iteration is set up so that p (l) and p^(L) approach the 

same distribution. 

Thus, for numerical purposes a sufficient number of terms must be retained in 
the series expressions for p^, and p a to adequately approximate each of 
these infinite series. Note that the numbers of terms required to approximate 
the series expressions for 9p a /3x and dp^/ax are likely to be different from 
the number required for p and p expressions, because the gradients vary more 
rapidly with position than the pressures. 



One further factor is worth noting. The series expressions for Bp^/dx and 
Bp^/Bx (Eq. B-l and B-2) suggest that fewer terms are required for positions 

on either side of the interface. The coefficients are reduced by the ratio 

_ ... . m7r(x-L)/w . _ , .. qrr(L+.-x)/w 

of sine terms which approaches e for x< L and large m or e - 11 

for x > L and large q. 
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APPENDIX C 


TWO-DIMENSIONAL BAFFLED CHAMBERS 


This appendix details the equations, logic and computer programs which were 
developed to described two-dimensional baffled thrust chambers with non-zero 
end-wall admittance boundary conditions. The rigid wall case is obtained by 
simply equating these admittance values to zero. 

Also included in this Appendix is a listing of the main computer program 
and the auxiliary programs used for generating pressure and velocity distri- 
butions and input files. Finally, a sample calculation is included for a 
rigid end-wall chamber containing a single baffle. 

To simplify the computer operations, the mathematical notation used here 
and in Appendix D is different from that used in the main text. A separate 
table of nomenclature for these appendices is included at the end of Appendix 
D, which also shows the relationships between the nomenclature used in the 
text and in the appendices. 


GENERAL EQUATIONS 


The general characteristic equation derived from the variational method for 
the two-dimensional baffled chamber shown in Fig. 1 is 



0 


at x = x 


(C-l) 


C-l 



It may be observed that this is simply an equation for conservation of the 
energy flux across the interface between the main chamber and the baffle 
compartments for the rigid-wall case. Corresponding to any pressure 
gradient at the interface,the pressures in the different regions are given 
by the integral equation corresponding to the wave equation,, i.e.. 


(x,y) = 


P' G ' (x,y) = + 


/ M+1 (U) at>( M ) 

I ' (x, y|x 1 , y 1 ) (x 1 ,y‘ ) dy> 

■'ty 

f Vp- (Q) _(^) 

I G ' (x, yjxp y') — (Xpy') dy' 

J Yo 


(C-2a) 
( C -2b ) 


( ii ) ( o ) -Kh 

where G ' and G are the Green's functions in the g baffle compartment 

and the main chamber respectively. The pressures obtained from Eq. C-2a 

and C-2b are not equal at the interface unless the pressure gradient at the 

interface used for the integration corresponds to the exact solution of the 

wave equation and the appropriate boundary conditions. 


Because the exact solution is unknown at the outset, an initial form for this 
pressure gradient is assumed and an iteration scheme is used to improve the 
accuracy of the initial form. This iteration scheme proceeds as follows. 

First, some form is assumed for the pressure gradient at the interface: 

3x (C-3) 

Then the above integral equations are used to obtain (x,y) and 

p' ' ; (x,y) directly. Since, for the exact solution, the pressures must 

match identically at the interface (x = x-^)> a new approximation for the 
baffle compartment pressure profile at the interface is obtained by equating 

the pressures ; 
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(C.-b) 


p (m)(C+ 1 )(^y) = p ( °) ) ( x^jjr) 

A new pressure gradient is then obtained by differentiating the corresponding 
integral equation and the iteration sequence is complete. The iteration may 
be repeated as many times as desired. 


Assumptions and Definitions 

The assumptions and definitions used to obtain the above algebraic forms of 
the characteristic and iteration equations are presented here for reference. 


The approximation to the pressure profile on either side of the interface 

is assumed to be represented by an eigenfunction expansion, i.e.. 
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where the orthonormal eigenfunctions are 
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the corresponding eigenvalues are* 
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The Kronecker delta is defined by: 

5 - f 1 ir. = n 

m^n = 

SO m / n 

Integrals of products of the eigenfunctions are then: 
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The Green's function for each region has been represented by an eigenfunction 
expansion also, 
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Final Equations 


When series expressions are used to represent the interface pressure profiles 
in both the main chamber and the baffle compartments, the iteration scheme 
proceeds as follows : 
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With the same series expression representations, the general characteristic 
equation reduces to an algebraic equation which can be solved numerically, 

i ,e . , 
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COMPUTER PROGRAMS 


The computer programs used to solve for the fundamental acoustic modes of 
the two-dimensional baffled chamber with gain/loss end -wall boundary condi- 
tions are described in this section. The computer logic is illustrated by 
the flow chart shown in Fig. 13 . All programs are written in 0E-400 Series 
Fortran* using complex arithmetic. The individual programs, functions and 
relation to equations and lines in the program are detailed below. 


B2D444 - Main Program (See Flowchart C-l) 

This program is used to obtain the frequency and damping coefficient (eigen- 
value) as well as the series expansion coefficients for the interface pres- 
sures and axial velocities for a particular two-dimensional baffled chamber 
with gain/loss end-wall boundary conditions. Complex quantities are noted. 

A full sample calculation is given at the end of this Appendix. 


a. Input 1: Line 1240 

MUBAR = b = Wo. of baffle compartments 
YO « y Q 
YMQBAR= yp 


b. Perform y Integration: Lines 1340-1760 : Eq. C-ll 

FISTIM(MOD,MTOD,HYnD)' - r " 


c . Input 2 : Line I 790 

XO * Xq 

XI = Xq 

X2 = x£ 

BETAOP = 0 (Complex) 

BETA2L = (3 (Complex) 

x 2- 


*CPB-1473, GE-400 Series Time-Sharing FORTRAN manual. General Electric Company, 
Information Systems Division, 13430 North Black Canyon Highway, Phoenix, 
Arizona, 1970. 
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Figure 13. Flow Chart for Program to Analyze Two-Dimensional Baffled Chambers 
The Accompanying Description of the Computer Programs is Keyed to 
This Flowchart Through the Letters a...k. “ 























d. Input 4: Line 1830 

MYOVMC = No. of antisymmetric terms retained in main chamber series 
NYUVMC = No. terms retained in baffle compartment series 
ITT - No. iterations performed 


Input 5 : Line 1880 

I TYPE = 0 for step calculation of characteristic function 

1 for root finding calculation of characteristic function 

Input 6: Line 1910: ( TTYPE - 0) 

EIYDST = k(y^--y 0 ) Initial guess (Complex) 

EIYDSP = k(y£-y 0 ) Increment (Complex) 

NEIYDS = No. of calculations 

IDETAL = jo for simple printout of characteristic function 
jl for detailed printout of series expansion 
l coefficients and pressure/axial velocity profiles 


Input: Line 1050: ( TTYPE = 0)( IDETAL = l) 


IPBSUM 


IPASOM 


IPVBUM 


IPVAOM 


l 0 for no print 

[ 1 for printout 


BSUM(MUD,NYUD) = 


asom(myod) 


= a ( o )(0 


VBUM(MUD,NYUD) = (x. ) 

U V Hlly X 

VAOM(MYOD) = (x, ) 


Input: Line No.: ( ITYPE =0)( IDETAL = l) 

MPBSUM ] ! 


MPASQM 


MPVBUM 


MPVAOM 


— 0 for no print 

1 for print 


(M)(c), , 

P v y) 

P (0)(C) (x 1 ,y) 

9p(^ )(* )(x 1 ,y)/3x 

9p(°)(^ )(x L ,y)/dx 


f. Input 7: Line 2000: ( ITYPE = l) 

EIYDST = k(yp -y Q ) Initial guess (Complex) 

EIYDFD = k(yjy-y o ) Finite difference (Complex) 

NEIYDM = No. of Calculations 

CHFNTL = Tolerance on characteristic function final value 
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Definition: Lines 2150-2730 ‘ Eqns. C-7a>t>; C-l4a,b; 

C-15a,b; Cl6a,b 


JT 

= 0 


ayom(myod) 

= (0) 

- a .. _ 


AX0V 

eh 
^3 * 

111 

(Complex) 

DELLOV 

= 6(0) 

x,my 

(Complex) 

hlom(myod) 

111 

^ 0 

>r 

H 

(Complex) 

bsijm(mud,nyud) 

III 

o' 

0 

( Complex) 

ayum(nyud) 

= a 

Y,ny 

( Complex) 

AXUV 

. «■ M 

x.n 

' y 


DELRUV 

hi 

Sotr 

( Complex) 

hrum(nyud) 

1 — I 
hi 

( Complex) 


Calculate Iteration Dependent Variables : 

Sqns. C-17, 0-18 

vbum(mud,nyud) 5 b^^ (x 1 ) 


blom(myod) 

asom(myod) 

vaom(myod) 


“y 

. B (0)(«) 

"V 

s a (o)(«) 

■ 4 ° )(£) 


Lines 2760-2950; 


( Complex) 
( Complex ) 

( Complex) 
( Complex ) 


Perform Iteration: /Lines 2960 - 3 HO: Eqn. C-l8e 
JI = JI +1 

BSUM(MUD,WYUD) = b ( M )( c+1 ) 

T 
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1. Calculate Characteristic Function: Lines 3140-3320: Eq. C-19 


SUM1 


SUM? 


Y ('|a(°K «)| 2 s ( 0 ) lx 

4^1 -v 1 \ ( Vj 


I 

H-l 


( Complex) 


2 (K )w Ug<*o) < 

V 1 


Comnlex) 


CHARFN = SUM1 - SUM2 


( Complex.) 


m. General Printout: Line 3360 
eignyp = k( y _-y o ) 

SUM1 

CHARFN 

SUM2 


( Complex) 

( Complex ) 
( Complex) 
( Complex) 


n. Detailed Printout: Lines 339^ - 36IO 

(According to Inputs g,h) 

Calls PV2DBC , PV2DMC for pressure/axial velocity profiles 


o. Increment Eigenvalue: Lines 3630 - 3950 

ITYPE = 6 - EIGNYD - EIGNYD + E1YDSP (Complex) 

ITYPE = 1 - EIGNYD = EIGNYD - CHAREN (Complex) 

‘ d( CHARFN) 
d( EIGNYD) 


p. Return to Options: Lines 3970-4010 

Input 1 IG0TO = 1 
"2 2 
" 4 4 

" s 5 

" 6 6 

" 7 7 
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Sub to u t ines 


FV2DBC - Baffle Compartment Interface Pressure/Axial Velocity Profiles 

■bis program takes coefficients of pressure/axial velocity series expressions 

and generates the appropriate profiles at the interface on the baffle 

c ompa rtment s i de . 

a. Input: Line 5000 (Variable List) 

NYPWT = Wo. of equally spaced y points for printout 
IMUD = Compartment at which printout should stop 
(NYPNT = 0 for return to main program) 

b. Printout: Line 5^20: Eqn. C-5a 

Value of pressure or velocity at y location 
y Location 
Compartment number 

: V2DM0 - Main Chamber Interface Pres sure/ Axial Velocity Profiles 

This program takes coefficients of pressure velocity series expressions and 

generates the appropriate profiles at the interface on the main chamber side. 

a. Input: Line 6000 (Variable List) 

WYPWT = Wo. of equally spaced y points for printout 
YST/p = y position at which printout should stop 
(NYPWT = 0 for return to main program) 

b. Printout: Line 632 O: Eq. C-5b 

Value of pressure or velocity at y location 
y location 
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B2D444 


0100C TWO-DIMENSIONAL BAFFLED CHAMBER MAIN PROGRAM 
1GUO DIMENSION FINTTMC 4, 24, 24) 

1010& , AY OM ( 2 4 ) >AY UM (24) 

1020& , HLOMC 24) ,HRUMC 24) 

1030& , AS0MC24) , VAOM ( 24) , BSUM C 4, 2 4) , VBUM ( 4, 2 4) , BLOW (24) 

10 40 COMMON MUBAR, YD, YDU, AYOM,AYUM,PI 

10 SO COMPLEX BETAOP , BETA2L, El YDST, E I Y DSP , El YDFD, El GNYD, EIYDQR, F.YDQR2 
1060& , AXOV, ARGO, RTARGO, DELLOV, CTANO, AXUV, ARGU, RTARGU, DELRUV, ctanu 
1 0 70& , HLQM, HLOV, HRUM,HRUV 

1080& , ASOM, ASOV, ALUM, ALUV, VAOM, BSUM, BSUV, BLOM, BLOV, VBUM 
1090& , TM, SERSUM, SUM1 , SUM2, CHARFN, CHARFO, CHARFL, CHARFU 
1100& , CHARFL, CHARFU, RTARGO, RTARGU 

1110 11 FORMAT ( / /''T W0 - DIMENSI 0NAL BAFFLED CHAMBER") 

1120 12 F0RMATC //T2, "MUBAR='M2/ 

1130& T4, "Y0=",F5.2,T1 7, "YMUBAR=", F 5 . 2, T34, "X0=", F S . 2, T47, "X 1 = ", F 5 . ?, 
1140& T60 , "X2="jF 5*2/ 

1 1 50& Til, "8ETA0P=",2F5.2, T31 , "BETA2L=", 2F5.2) 

1160 14 F0RMATC/T2, "TRANSVERSE M0DE NUMBER", 13) 

1170 16 F0RMAT(/T2, "MYOVMC=", I2/T2, "NYUVMC=", I2//T2, "JI T=", 12) 

1180 18 F0RMATCT16, "El GNY D", T52, "SUM 1"/T 1 6, "CHARFN", T52, "SUM2") 

1190 19 F0RMATCT1O, 1PE16. 10,",", 1PE1 7. 10) 

1200 20 F0RMAT(2( 1PE16.10,",", 1PE16.10," : ", 1 P E 1 6 . 1 0, ", ", 1 P E 1 6 . 1 0 / ) ) 
1210 PI=3. 1 41 5926536 

1220C MUBAR IS THE NUMBER 0F COMPARTMENTS 
1230.1001 PRINT," 1 I NPUT: MUBAR, YO, YMUBAR" 

1240 INPUT, MUBAR, YO, YMUBAR 
12S0 IFCMUBAR.LT. 5) G0 T0 10015 

1260 PRINT, t, "MUBAR MUST BE LESS THAN 5--TRY AGAIN" 

1270 G0 T0 1001 
1280 10015 CONTINUE 
1290 FMUBAR=FL0AT (MUBAR) 

1300 YD=YMUBAR-YO 
1310 YD2=YD**2 
1320 Y DU=YD/F MUBAR 
1330C 

1340C Y INTEGRAL CALCULATION 
1350C 

1360 MY 0 1 NT = 2 4 

1370 NYUI NT=24 

1380 DO 910 MUD=1, MUBAR 

1390 F I NTTM (MUD, 1 , 1 )= 1 • O/SQRT (FMUBAR) 

1400 DO 901 NYUD=2, NYUINT 
1410 F I NTTM (MUD, 1,NYUD)=0.0 
1420 901 CONTINUE 
1430 FMUD=FLOAT (MUD) 

1440 YUL0=YDU*CFMUD-1 .0) 

1450 YUP 1 L0=YDU*FMUD 
1460 DO 904 MY0D=2, MY"OI NT 
1470 FMYOD=FL0AT(MYOD) 

1480 FMY0D=2 .0*CFMY0D- 1 .0) 
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B2 D - / j A A CONTINUED 


1490 AL Y 0 V = ( F M Y 0 D - 1 .0)*PI/YD 

IbOO Ah:GUU=ALYOV*YULO 

1 b 1 0 A R G UP 1 = AL Y QV*YUP 1L0 

1 520 0 = ( SORT C 2 .Q*FMUBAR) ) /< ALYOV*YD) 

1530 fv AL UE= C* C 5 1 N ( ARGUP 1 ) -SINCARGUU) ) 
i 540 IFC ( ABS(FVALUE) ) -LT. ( 1 -0E-6) ) FVALUE=0.0 
1 b 5 0 F I N T T M ( MUD* MYO D> 1 ) = F V A L IJ E 
1 560 u£ 903 NY UD = 2> NYUINT 
1 o 7 j F N Y U D = F L GAT (NY U 0) 

I b fc 3 A L Y U V = ( F NY LJ D- 1 • 0 ) *P I / YDU 
1 b 9 J A L Y D I F = A B S C AL Y 0 V - AL Y U V ) 

ISDU IF < ALYDIF *LT • ( 1 «0E-6) ) G0 T0 902 

1610 G*CALY0V/<ALY0V**2-ALYUV**2> ) *SQRT( FMUBAR) *2 • O/YD 

1 6 20 A C = ( F NY OU- 1 • 0 ) *P I 

1630 F VALUE = C* ( COS ( AC ) *SI N( AR GUP 1 ) - S I N ( AR GUU) ) 

1640 IF ( ( ABSCFVALUE) ) «LT. C 1 .OE-6) ) FVALU£=0.0 
1 6 5 0 F INI T i v i ( MUD, MYQD, NYUD) =FVALUE 

*660 GO i 0 903 

1 6 7 J 9 0 2 CONTI N IJ E 

1 6 8 0 C=1 • 0 /SQRTC FMUBAR) 

1 6 9 0 FVALUE=C*C0SCARGUU> 

1700 IFC ( ABSCFVALUE) ) .LT. ( 1 .OE-6) ) FVALUE=0.0 

1 7 1 0 F I N T T Fi ( M UD, MYO D, NY UD ) = F VALUE 

1 720 903 CONTINUE 

1 - 3 0 9 0 4 C 0 N T I N U E 
1740 910 CONTINUE 
1 1 b 0 C 

1 7 6 0 C END Y INTEGRAL CALCULATION 

1 7 70 C 

1780 1002 PRINT," 2 INPUT:X0,X1,X2,BETA0P, BETA2L" 

1790 I NPUT, XO, X 1 , X2, BETAOP, BETA2L 

IBOO 1003 CONTINUE 
1 B I 0 Nib 0 D E = 1 

1820 1004 PRINT," 4 I NPUT sMYOVMC, NYUVMC, JI T" 

1 8 3 0 I N ? U T , MYOVMC, NYUVMC, JI T 

IB 40 IFC (MYOVMC -LE.MYOINT) -AND. (NYUVMC. LE. NYUINT) ) 60 T0 1005 

1850 PRINT, t, "MYOVMC & NYUVMC MUST BE LESS THAN 25--TRY AGAIN" 

1 8 6 C G 0 T £ 10 0 4 

1870 1005 PRINT," 5 I NPUT: I TYPE" 

1 8 8 C INPUT, I TYPE 

i o j 0 i U tL I ML. — U 

1900 I F ( I TYPE) 9999, 1006, 1007 

1910 1006 PRINT," 6 I NPUT : El YDST, El YDSP , NEI YDS, I DETAL" 

1920 I NPUT, El YDST, El YDSP, NEI YDS, I DETAL 

1930 IFC I DETAL .NE. 1 ) G0 T0 1008 

1940 PRINT, " INPUT: IPBSUM, IPASOM, IPVBUM, IPVAOM" 

1950 INPUT, IPBSUM, IPASOM, IPVBUM, IPVAOM 

1960 PRINT,” I NPUT : MPBSUM, MPASOM, MPVBUM, MPVAOM" 

1 9 70 I NP U T , MP B S UM , MP A SOM , MP VB UM , MP VAO M 

1 9 8 0 G 0 T 0 1008 


c-l6 



B2D444 C0NTINUED 


1990 1007 PRINT^" 7 INPUT: EIYDST* El YDFD> NEI YDM, CHFNTL" 
2000 INPUT-* El YDST., El YDFD* NEIYDM* CHFNTL 
2010 1008 CONTINUE 
2020C 

2030C PRINT INPUT 

20 40 C 

2080 PRINT 11 

2060 PRINT 12>MUBAR*Y0jYMUBAR*X0»X1*X2»BETA0P>BETA2L 

2070 PRINT l At NMO DE 

2080 PRINT 16,MY0VMC,NYUVMC* JIT 

2090 PRINT 18 

2100 NEI Y DD=0 

21 10 NCNTRL=- 2 

2120 EIGNYD=EIYDST 

2130 1 CONTINUE 

21 40 C 

2150C DEFINITIONS 

2160C 

2170 JI=0 

2180 EIYD0R=EI GNYD/YD 
2190 EYD0R2=EIYD0R**2 
2200 AY OM Cl) =0.0 
2210 AY UN C 1 )=0.0 
2220 DO 30 MY0D=2 jMY0VMC 

2230 FMYO D=FL0ATC MYOD ) 

2240 FMY0D=2.0*(FMY0D-1 .0) 

2250 AYOV=CFMYOD- 1 *0)*PI/YD 
2260 AYOMCMYOD)=AYOV 
2270 AY0V2=AY0V**2 
2280 RTARG0=EYD0R2- AY0V2 
2290 AXO V=CSQRT ( RTARGO ) 

2300 ARG0=AX0V*CX1-X0> 

2310 DELLO V= C 0 • 0* 1 «0)*EIYD0R*BETA0P/AX0V 
2320 CTANO=CS INC ARGO) /CC0SCARGO) 

2330 HLOV= - AXO V* ( CTANO+ DELLO V) /C 1 .0- DELLOV*CTANO ) 

2340 HLOMCMYOD) =HLOV 

2350 30 CONTINUE 

2360 HLOM Cl) =0 * 0 

2370 D0 32 MUD= 1 > MUBAR 

2380 D0 31 NYUD= 1 » NYUVMC 

2390 BSUMCMUD* NYUD) =0 -0 

2400 31 CONTINUE 

2410 32 CONTINUE 

2420 IF C MUBAR . NE * 1 ) G0 T0 33 

2430 BSUMC 1,2)=1 .0 

2440 G0 TO 39 

2450 33 I F C MUBAR • NE • 2 ) G0 T0 34 
2460 BSUMC 1 , 1 ) = +l .0 
2470 BSUMC 2* 1 ) = - 1 - 0 
2480 G0 TO 39 
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B2D444 CONTINUED 


2490 34 I F ( MUBAR • NE * 3 ) G0 T0 35 

250 0 BSUMC 1 > 3 )=+l .0 
23 1 0 BSUMC 2* 1 )= 0.0 
2320 BSUMC3* 1 )=- 1 *0 
2 b 3 0 G 0 T 0 39 

2340 33 IF C MUBAR • NE • 4) GG T0 36 

25 50 B S UM C1*1)=+1»Q 
2360 B S UM ( 2 > 1 ) = + 0 • 5 
2370 BSUM ( 3* 1 ) =-0 • 5 
252 0 BSUMC 4* 1 )=- 1 .0 
2590 36 CONTINUE 
2600 39 CONTINUE 

26 10 00 40 NY IJD 3 1 » NYUVMC 
26 2 0 F NYUD=FL0ATCNYUD) 

2630 AY UV = C F NYUD- 1 .0)*PI/YDU 
2640 AYUMCNYUD) =AYUV 
2650 AYUV2=AYUV**2 

26 6 0 R T A R G U = E Y DO R 2 - AY U V 2 
2670 AXUV'-CSQRTCRTARGU) 

2680 ARGU=AXUV*CX2-X1 ) 

2690 DELRUV=(0 -Oj 1 .0) *EIYD0R*BETA2L/AXUV 
2 7 00 C T A N U = C S I N C A R G U ) / C C 0 S C A R G I J ) 

2710 HRUV = AXUV* C CTANU+DELRUV) /< 1 .0- DELRUV*CTANU> 

2 7 20 H R UM C N Y U D ) = HR U V 
2730 40 CONTINUE 
2740 45 CONTINUE 

27 500 

2 7 6 0 C BLOW CALCULATION 
2 7 7 0 C 

2780 DO 90 MY0D=2> MYOVMC 

279 0 S El R SUM =0.0 

2800 D0 80 MUD= 1 » MUBAR 

2b 10 D0 60 NYUD=1> NYUVMC 

2820 TM = BSUM C MUD j> NYUD) *HRUM C NYUD) 

28 3 0 V B UM (MUD.* NYUD) =TM 

28 40 TM= TM*F I NTTM C M UD* M YO D* NY UD ) 

28 50 SERSUM= SERSUM + TM 
2860 60 CONTINUE 
2870 80 CONTINUE 

2880 BLOM ( MYOD) =SERSUM 
2890 VAOM C MYOD) =SERSUM 

29 00 ASOM C MYOD) =SERSUM/HLOM CMYOD) 

2910 90 CONTINUE 

29 20 3L0MC 1 )=0 .0 

29 30 VAOM C 1 ) =0 • 0 

29 40 ASOM Cl) =0 « 0 

2950 IFC JI .EQ • JIT) G0 T0 150 

2960 JI=JI+t 

29 70 C 

2980C ITERATION EQUATION CALCULATION 
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B2D444 CONTINUED 


2990C 

3000 D0 140 MUD= 1 * MUBAR 
3010 D0 120 NYUD= 1 * NYUVMC 
3020 SERSUM=0 • 0 
3030 D0 110 MY0D=2*MY0VMC 
3040 TM=AS0M C MYOD) 

3050 TM=TM*FINTTMC MUD* MYOD* NYUD) 

3060 SERSUM=SERSUM+TM 
3070 110 CONTINUE 

30 8 0 BSUM CMUD*NYUD)=SERSUM 
3090 120 CONTINUE 

3100 140 CONTINUE 
3110 GO TO 45 
3120 150 CONTINUE 

31 30C 

3140C CHARACTERISTIC EQUATION CALCULATION 
31 50C 

3160 SERSUM = Q • 0 

3170 DO 170 MY0D = 2*MY0 VM C 

3180 ASOV=ASOMCMYOD) 

3190 TM= C AS0V**2) *HL0MCMY0D) 

3200 SERSUM=SERSUM+TM 

3210 170 CONTINUE 

3220 SUM 1 =SERSUM 

3230 SERSUM=0 • 0 

3240 DO 200 MUD= 1 * MUBAR 

3250 DO 190 NYUD= 1 * NYUVMC 

3260 BSUV=BSUMCMUD* NYUD) 

3270 TM=(BSUV**2)*HRUMCNYUD) 

3280 SERSUM=SERSUM+TM 
3290 190 CONTINUE 
3300 200 CONTINUE 
3310 SUM2=SERSUM 
3320 CHARFN = SUM 1 - SUM 2 
3330C 

3340C GENERAL PRINTOUT 
3350C 

3360 IFCNCNTRL. EQ .-2) PRINT 20* El GNYD* SUM1 » CHARFN* SUM2 
3370 IF ( I DETAL • NE • 1 ) GO TO 260 
338 OC 

339 OC DETAILED PRI NT0UT- - BSUM* ASOM* VBUM , VAOM 
3400C 

3410 IFC C IPBSUM+MPBSUM) .NE.O) PRINT* t2* ”BSUM”, t 
3420 IFC IP BSUM • NE • 1 ) GO TO 241 

3430 PRINT 1 9 * ( C BSUMC MD* NY ) * MD= 1 * MUBAR) * NY= 1 * NY UVMC ) 
3440 241 CONTINUE 

3450 IFC MPBSUM • EQ • 1 ) CALL PV2DBCCBSUM* NYUVMC) 

3460 IFCCIPASOM+MPASOM) .NE.O) PRI NT* *2* "ASOM", t 
3470 IFCIPASOM.NE'.l) GO TO 243 
3480 PRINT 19* CAS0MCMY)*MY=1* MYOVMC) 
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3490 243 CONTINUE 

3b0u I F ( Op ASOOi • EQ • 1 ) CALL PV2DMCCAS0M*MY0VMC> 

351 C IF C C IPVBUM+MPOBUM) .NE.O) PRINT* tQ, "VBUM"* t 

3320 IF( IPVEUM.NE. 1 ) G0 T0 245 

3530 PRINT 19* ( < VBUtMCMD* NY) * MD= 1 * MUBAR) * NY = 1 * NYUVMC) 

3 3 4 > 4 3 , 4 N T I N I J L 

3330 IF (MPVBUM.EQ. 1 ) CALL PV2DBCC VBUM* NYUVMC) 

3560 I F (IIP OAQM+MP VAOM) • NE • 0 ) PRINT* 1 2* " VAOM" » t 

33 70. I F C I p y A. 0 M « N E * 1 ) G0 T0 247 

33 3 u PRINT 1 9 > CVAOMCMY)*MY=1*MYOVMC) 

35 9 i 247 CONTINUE 

3600 I F ( OF 97000! .£0.1) CALL PV2DMCC VAOM* MYO VMC ) 

36 10 26(1 CONTINUE 
3620 C 

36 3 0 C CHANGE EIGNYD--STEP CHANGE 0R R00T FINDER 

36 3)0 0 

36 3 0 N E 1YE||=NEIYDD+1 

36 60 IF ( I TYPE « EQ • 1 ) G0 T0 300 
3670 FI GNY D=EI GNYD+EI YDSP 

3660 I F (NEIYDD-LT.NEIYDS) G0 T0 1 

3690 GO 10 1009 

3709.7 300 CONTINUE 

3710 EIYDFR=EIYDFD 

3720 E I Y DF I =■>■ (0*0* 1 .0)*EIYDFD 

3730 IFCABSCEIYDFR) .GT. C 1 .OE- 10) > G0 T0 301 

3740 CHARF 1 =~ CO .0-* 1 .0)*CHARFN 

3733 CHARF 0 = + C 0 • 0 * 1 .0>*CHARFI 

3 760 GO T 0 30 2 

37 7(7 30 1 CONTINUE 

3780 IFCABSCEIYDFI) .GT.C1 .0E-10)) G0 T0 302 

3790 CHAFQ* R = CHARF N 

3800 C H A R F N = C H A R F R 

3810 302 CONTINUE 

36 2:7 NCNTRL = NC;NTRL+1 

3830 IF < NCNTRL) 310* 320*330 

3840 310 CHARFQ= CHARF N 

38 30 I: 7 C CAB3CCHARF0) .LT.CHFNTL) .0R. CNEIYDD. GT.NEIYDM) ) G0 T0 1009 

3860 £1 GNYD=EIGNYD-EIYDFD/2.0 

38 70 COO 20 | 

3830 320 CFARFL=CHARFN 
3890 FI GNYD=EI GNYD+EI YDFD 
3900 GO TO I 
3910 330 CHARFU=CHARFN 

39 20 E i GNY D = El GNYD- El YDFD/2 • 0-CHARF0*EI Y DFD/ ( CHARFU-CHARFL ) 

3930 NCNTRL= - 2 
39 40 GO T0 I 
3950 1009 PRINT* t2*t2 
39 6 0 C 

39 70 C WHERE T0 G0? 

39 8 0 C 


C-20 
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B2D444 CONTINUED 


3990 PRINT*" INPUT: IG0T0" 

40 0 0 INPUT* IG0T0 

4010 G0 T0 C 1001 * 1002* 1 003* 1004* 1005* 1006* 100 7* 1 008*9999) * I G0T0 
4020 9999 STOP 
4030 END 


C-21 
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PV2DBC 


50 0 0 SUBR0UT I NE P V2DBCC C0EF* NYUVMC) 

50 1 0 DIMENSION C0EF( 4* 24) > AYOM ( 24) * AYUM C 24) * EPSM C 24) * PVVAL C 50 ) * YVALC 50) 

5Q2C& j. MU DM (50) 

5030 C0MM0N MUBAR* Y D* YDU* AYOM* AYUM* P I 

30 40 COMPLEX C0EF 

5030 1 5 F0RMATC/5X* "Y- VARI ATI 0N") 

5060 10 F0RMATC2X* 1P2E1 7. 10* 110) 

50 70 Y DU02 = Y DU/2 • 0 

3080 DO 100 NYUD= 2* NYUVMC 

50 3 0 £ P S M < N Y U D ) = S Q R T < Y D U 0 2 ) 

3100 100 CONTINUE 
3110 EPSMC I ) = SQRTCYDU) 

3120 1 CONTINUE 

3130 PRINT* t2>" I NPUT : NYPNT* IMUD"* t* 

51 40 INPUT* NYPNT* IMUD 

3150 IF (NYPNT. EQ *0) RETURN 

3160 Y INC=Y D/NYPNT 

5170 NYPNT 1 =NYPNT+1 

3180 YDUS VE=Y DU 

5190 YDU=YDU+1 .OE-6 

5200 Y P 0 S = 0 « 0 

5210 DB 500 NY D=l* NYPNT 1 

5220 IF ( ( YP0S * GE • ( 0 • 0 ) ) • AND • ( YP0S *LE • YDU) ) MUD=1 
3230 1 F ( (YP0S.GT.YDU) .AND- (YP0S.LE* (2.0*YDU) ) ) MLID=2 
5240 IF( ( YP0S -GT. (2.0+YDU) ) . AND . ( YP 0S * LE . C 3 . 0*YDU) ) ) MUD=3 
5250 IF ( (YP0S .GT. (3 .0*YDU) ) .AND.C YP0S.LE. ( 4.0*YDU) ) ) MUD=4 

5260 MUDM ( NYD) =MUD 
32 7 0 F M UD = FL0AT (MUD) 

5280 SUM = 0 » 0 

5290 D0 400 NYUD=1 * NYUVMC 

5300 TM=(C0S(AYUM(NYUD)*(YP0S-(FMUD-1 .0)*YDU) ) )/EPSM(NYUD) 

5310 5UM=SUM+TM*C0£FCMUD* NYUD) 

5320 400 CONTINUE 
33 3 0 P V V ALC NYD) = SUM 
5340 YVALC NYD) =YP0S 
5350 YP0S=YP0S+YINC 
5360 N Y MAX 3 NY D 

5370 IF (MUD. EQ • IMUD) G0 T0 600 

5380 500 CONTINUE 

5390 600 CONTINUE 

5400 YDU-YDUSVE 

5410 PRINT 15 

5420 PRINT 10* (P VVALC NYD) * YVALC NYD) *MUDM C NYD) * NYD= 1 * NYMAX ) 

5430 G0 T0 1 
5440 END 
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1 


PV2DMC 


6000 
6010 
6020 
6030 
60 40 
60 60 
6060 

60 70 
6080 
6090 
6100 

61 10 
6120 
6130 
61 40 
61 50 
6160 

61 70 
6180 
6190 
6200 
6210 
6220 
6230 
6240 

62 50 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 


SUBROUTINE PV2DMC< C0EF>MYOVMC> 

DIMENSION C0EFC24) , AYOM < 24) , AYUM ( 24) , EP SM C 2 4 ) , P V VAL C 50 
COMMON MUBARjYD* YDU* AYOM, AYUM, PI 
COMPLEX CO EF 

1 5 FORMATC/SX^'Y-VARIATION") 

10 F0RMATC 2X> 1P2E1 7. 10) 

YD02=YD/2 .0 

DO 100 MYO D= 1 j MYOVMC 

EPSMCMYOD) =SQRTCYD02> 

100 CONTINUE 
EPSMC 1 ) =SQRT( YD) 

1 CONTINUE 

PRINTj *2j" INPUT: NYPNTj YST0P”# t* 

INPUTS NYPNT* YSTOP 

I F < NYPNT • EQ • 0 ) RETURN 

YINC=YD/NYPNT 

NYPNT1 =NYPNT+ 1 

YP0S = O .0 

DO 500 NYD= 1 .» NYPNT1 
SUM=0 • 0 

DO 400 MYOD= 1 » MYOVMC 

TM= ( COS CAYOMCMYOD)*YP0S>> /EPSMCMYOD) 

SUM= SUM+TM*C0EF ( MYOD) 

400 CONTINUE 
PVVALCNYD) =SUM 
YVAL C NYD) =YP OS 
YP0S=YPOS+YINC 
NYMAX= NY D 

IFCYP0S.GT.YST0P) GO TO 600 
500 CONTINUE 
600 CONTINUE 
PRINT 15 

PRINT 10* (PVVALCNYD)>YVAL(NYD)>NYD=1>NYMAX) 

GO TO 1 
END 


YVALC 50) 



B 2 0 4 4 4 


20 s 49 NAR JAN. 12>1971 


1 I NP UT s MUBARj Y0> YMUBAR 
? 2. -0.0 >20.0 

2 I NPUT : XO > X 1 > X2> BETAOP> BETA2L 
?Q «0» 2! 2 • 0 > 30 «0> 0 » 0 > 0 « 0 > 0 • 0 > 0 *0 

4 I NPUT : MYO VMC* NYUVMC-. JI T 
? 9 > 1 1 1 9 

5 INPUT: I TYPE 
?0 

6 I NP LI T ; E I Y DST > E I Y DSP > NE I Y DS * I DETAL 

? 2 . 10 * 0 . 00 * 0 . 10 > 0 . 00 > 5>0 



TWO-DIMENSIONAL BAFFLED CHAMBER 


MU BAR- 2 

Y 0 = 0.00 YMUBAR=20 *00 X0= 0.00 Xl=22.00 X2 = 30.00 

BETAOP= 0.00 0.00 BETA2L= 0.00 0.00 

TRANSVERSE MODE NUMBER 1 

MYO v.'MC= 9 
NY U VMC = 1 1 


EIGNYD SUM 1 

CHARFN SUM2 

2 « 1 0 0 0 0 0 0 0 0 E + 0 0 > 0 .OOOOOOOOOE + OO : 6 . 538 I 8 420 5E-0 4> - 2 . 02 78 1 6548 E- 1 5 
1 .3321 1 S098E-03* - 7 . 9 339 74209 E- 1 5 : - 6 . 782966 779E-0 4, 5.9061 57662E- 1 5 

2 • 2 3000000061 + 00 .» 0 . OOOOOOOOOE+OO : 9 . 9 4009 7568 E- 0 4> 8 • 8 7342 7550E- 1 4 
1 - 50 78 50 29 5E-03> -4.960516071 E- 14 : -5. 1 38 40 538 1 E-0 4> 1 . 38339 4362E- 1 3 

2. 300000000E+00> 0 .OOOOOOOOOE+OO : 1 .622256384E-02, 6 . 630450630E- 1 2 
-9 • 60 6 727823E-0 4> ~8 • 5548633 78E- 1 3 : 1 . 7 1 8 323662E- 02, 7. 485936967E- 1 2 

2 . 400 0 0 OOQQE+OO, 0 .OOOOOOOOOE+OO : 1 . 1 73579350E+00, 4 . 920248 1 3SE- 1 0 
5.93679 1 1 42E-02* 4 . 88 1 3269 62E- 1 1 : 1 . 1 1 421 1 438E+00, 4 . 432 1 1 5441 E- 1 0 

2* b00000000E+00> 0 .OOOOOOOOOE + OO : 1 .0 33 1 72 5 1 6E+02, 4 . 1 9 73 52 1 5 4E-08 
2. 5261 1 4527E+G1, 1 . 1 9471 4700E-08 I 7 . 80 56 1 0635E+0 1 , 3 • 00263 7454E-08 


INPUTS I G0T0 
? 5 


C-24 




5 INPUT: ITYPE 
?1 


7 INPUT: EIYDSTjEIYDFD>NEIYDM> CHFNTL 
?2»40>0«00j 1 *0E-5j0«00j20/ 1 « OE - 6 


b. find root precisely 


TWO-DIMENSIONAL BAFFLED CHAMBER 


MUBAR= 2 

YO= 0-00 YMUBAR=20 *00 

BETAOP= 0.00 0.00 


X0 = 0.00 X 1 = 2 2 • 0 0 

BETA2L= 0.00 0.00 


TRANSVERSE MODE NUMBER 1 


MYOVMC= 9 
NY UVMC= 1 1 


X2= 30 . 00 


J I T= 9 

El GNYD 
CHARFN 

2. 400000000E+00, 0 .OOOOOOOOOE+OO 
5.936791 161E~02> 4 . 88 1 3269 72E- 1 1 


2.388201 239E+00, 
1 • 8335888 10E-02> 

2*380035851E+00j 

4.618831696E-03* 

2.376l93064E+00> 
6 .675558889E-04, 

2. 37542420 4.E + 00, 
2*22393664 4E-05> 

2«375396778E+00j 
2. 728484105E-08* 


0 .OOOOOOOOOE+OO : 
2.204684692E-1 1 : 

0 .OOOOOOOOOE+OO : 
1 . 191939794E-1 1 : 

0 .OOOOOOOOOE+OO : 
8 .65942891 4E-12 : 

0 .OOOOOOOOOE+OO : 
8 .097440466E-12 : 

D .OOOOOOOOOE + OO : 
8 .077901 597E-12 : 


1 . 1 7357935 1 E+OO* 
1.1 1421 1 440 E + 00 > 

6.999842057E-01 j 
6.8164831 76E-0 1 

4.901 332498E-0 1 * 
4.8551 44181 E-01 > 

4.1 4600360 1 £-0 1 > 
4*1 39328042E-01 , 

4.009582986E-01 > 
4.009360593E-01 ^ 

4*00480 1 1 20E-0 1 * 
4.004800847E-01 * 


SUM 1 
SUM 2 

4.920248 1 44F- 1 0 
4.4321 1 544 7 E - 1 0 

2*9 4444486 4E- 10 
2. 7239 7639 5E- 10 

2.0 660 40 71 BE- 10 
1 . 9 4 6 8 4 6 7 3 9 E - 1 0 

1 . 749232 59 4E- 1 0 
1 . 66263830 5E- 1 0 

1 . 69 19 69 56 8 !•> 1 0 
1 .61 09951 63E- 10 

1 .689962101 F> 10 

1 . 609 1 8308 5E- 1 0 


? 5 


INPUT: IG0T0 




6 I N P ! J T:EIYDSTj>EIY DSP » NEI YDS* I DETAL 

?2. 37539 67 78 >0*00^0. 00 >0. 00 > b 1 


INPUT: IPBSUM, IPASOM, IP VBUM, IPVAOM 

? 1 * 1 > 0 > 0 

I NP UT : MP 8SUM>MPAS0M>MPVBUM>MPVA0M 

7 1 » 1 j > 0 jO 


c. check pressure match 


TWO-DIMENSIONAL BAFFLED CHAMBER 


M U BA R = 2 

Y0= 0.00 YMUBAR=20.00 

BETAOP= 0.00 0.00 


X0= 0.00 X 1 =22.00 

BETA2L= 0.00 0*00 


X2=30.00 


TRANSVERSE MODE NUMBER 1 


MYOV, C= 9 
N Y U V M C = i 1 

JIT- 9 

E1GNYD SUM 1 

CHARFN SUM 2 

2 • 37S396 778E + 00 j> 0 .OOOOOOOOOE + OO : A .00 480 1 1 53E-0 ! * 1 . 6899 62 1 1 5E- 1 0 
2 » 7432 1 79 1 9E-08j> 8 . 0 7790 1 728E- 1 2 : A . 00 48008 79 E-0 1 ■> 1 • 609 1 83098E- 1 0 




BSUM 


1 .277951 188E+00, 
-l .277951 188E+00.* 
4.081 312083E-01 * 
4.081312083E-01> 
-1 . 1 5641 1 098E-0 1 > 
1 . 1 5641 1098E-01 * 
6 .061 904496E-02> 
6-061904497E-02, 
-4 .094246570E-02, 
4.094246571 E-02, 
3. 196590077E-02.> 
3. 1965900 78E-02* 
-2 • 782020370E-02, 
2 . 782020371 E-02^ 
2 .725257004E-02, 
2.725257005E-02* 
-3.578051283E-02, 
3.578051283E-02, 
2 .0449 1 6 720 E- 02 > 
2.04491 6720E-02, 
-1.491 370772E-02, 
1 .491 370772E-02, 


2.51 4920637E- 10 
-2*51 492Q6 37E - 1 0 
7 . 5835621 93E- 1 1 
7.583562194E-1 1 
-2.504733347E- 1 1 
2. 5047333 48 E- 1 1 
1 .3857671 24E- 1 1 
1 .3857671 24E- 1 1 
-9 .3448191 51 E- 12 
9 .3448191 53E- 12 
7 .077730648E- 1 2 
7.077730649E- 12 
-5. 85989258 7E-12 
5.859892588E-12 
5. 3 48 733180 E- 12 
5* 3 48 733182 E- 12 
-6.240244326E-12 
6.240244326E-12 
3* 727 40 8854 E- 12 
3. 7274088 54E- 1 2 
-2.759161578E-12 
2. 7 59 161578 E -12 


INPUT:NYPNT.» IMUD ? 40i 2 


Y- VARIATION 
. 5.442445646E-01 
5 . 493694995E-01 
5.537478926E-01 
5.447328838E-01 
5.300293430E-01 
5.213012463E-01 
5.1 42940501 E-01 
5.012466743E-01 
4.8446781 47E-01 
4.653907526E-01 
4.406395530E-01 
4.1 63369 593E-0 1 
4.0 10044394E-01 
3.813387165E-01 
3.388607584E-01 
2.91 563041 4E-01 
2. 727839578E-01 
2.635908060E-01 
2.064240570E-01 
1 .063659 120E-01 
5.372469061 E-02 
- 1 *0636572 43 E-01 


y 

O.OOOOOOOOOE+OO 
5.000000000E-01 
1 .OOOOOOOOOE+OO 
1 • 500000000E+00 
2. OOOOOOOOOE+OO 
2. 500000000E+00 
3. OOOOOOOOOE+OO 
3. 50 0000000 E+00 
4. OOOOOOOOOE+OO 
4. 500000000E+00 
5. OOOOOOOOOE+OO 
5.500000000E+00 
6. OOOOOOOOOE+OO 
6. 500000000E+00 
7. OOOOOOOOOE+OO 
7. 500000000E+00 
8. OOOOOOOOOE+OO 
8 • 500000000 E+00 
9. OOOOOOOOOE+OO 
9 . 500000000E+00 
1 .000000000E+01 
1 « 0 50000 00 OE+O 1 


compartment 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 


P ( " )( ' <, (x 1 ,y) 


INPUT: NYPNTj IMUD ?0,0 
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A 5 U ri 


0 .OOOOOOOOOE+OO, 
I *831 161819E+QO, 
-2 .940 8 986 42 E-01 , 
1 .2 379 9 79 06 £-0 1 , 
-6.96663361 IE-02, 
4.431 1 3 4790 E- 02, 
-3.072577357E-02, 
2.261 1269 12E-02, 
- 1 • 739 1 02532E-02, 


0 .OOOOOOOOOE+OO 
4. 100003583E- 10 
-5.471 580 740 E- 1 1 
2-459 738929 E- 1 1 
- 1 .46531 3026E- 1 1 
1 -001 758054E- 1 1 
-7.451021271 E- 1 2 
5.8823 8 0959E- 12 
-4. 872570233 E- 12 


Y- V/AkIATI ON 
5 • 0 9 7 6 2 7 7 2 6 E - 0 1 
5 * 10 79 49 39 6 E- 01 
:.>*! 1 4 54 59 1 0E-01 
5.0821 51 431 E-01 
4.0071951 33E-01 
4 * 9 2 9 3 3 8 2 6 0 E - 0 1 
4.87432031 5E-Q1 
4 . 8 1 6 0 0 5 2 9 2 E - 0 1 
4. 7 1 0032229 E-0 1 
4.543786531 E-01 
4 . 393472389E-01 
4 . 2 6 0 9 5 8 3 4 5 E - 0 1 
4. 121 3399 43E-01 
3 *908 1 3 400 6 E-01 
3*610 5 50 300 E~01 
3. 296 526929 E-01 
3.0 18264907E-0 1 
2 .69774709 1 E-01 
2* 1 4261 848 5E-01 
1 -21 5 4 39983 E ~01 
2.231341 1 44E- 1 1 
- 1 .21 54 3 9982 El - 01 


I NP UT : NYPNT, Y STOP 


y 

0 .OOOOOOOOOE+OO 
5. 0000 00000 E-01 
1 .OOOOOOOOOE+OO 
1 • 500000000E+00 
2. OOOOOOOOOE+OO 
2. 500000000E+00 
3. OOOOOOOOOE+OO 
3 *50 0000000 FI + 00 
4. OOOOOOOOOE+OO 
4. 500000000E+00 
5. OOOOOOOOOE+OO 
5. 500000000E+00 
6 .OOOOOOOOOE+OO 
6 . 500000000E+00 
7. OOOOOOOOOE+OO 
7. 500000000E+00 
8- OOOOOOOOOE+OO 
8 « 500000000 E + 00 
9 .OOOOOOOOOE+OO 
9 . 500000000E+00 
1 .OOOOOOOOOE+Ol 
1 .050000000E+01 


? 40 , 1 0 . 5 


p (0)(i)( Xl> y) 


INPUT : NYPNT, YST0P ?0,0. 


NPUTs IG0T0 


RUNNING TIMES 93.9 SECS 1/0 TIME : 22*5 SECS 
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APPENDIX D 


THREE DIMENSIONAL BAFFLED CHAMBERS 

This Appendix details the equations, logic and computer programs which were 
developed to describe three-dimensional baffled thrust chambers with non- 
zero end-wall admittance boundary conditions. The rigid wall case is obtained, 
by simply equating these admittance values to zero. 

Also included in this Appendix is a listing of the main computer program and 
the auxiliary programs used for generating pressure and velocity distribution 
programs used for generating input files. The necessary Bessel function sub- 
routines are also included. Finally, a sample calculation is included for a 
rigid end-wall chamber containing a three-bladed baffle. 

To simplify the computer operations, the mathematical notation used, hare and 
in Appendix C is different from that used in the main text. A separate table 
of nomenclature for Appendices C and D is Included at the end of this Appendix, 
which also shows the relationships between the nomenclature used in the text and 
in the appendices. 

This Appendix is written so that it may be used independently of Appendix C 
and is therefore somewhat redundant. 
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General Equations 


The general characteristic equation derived from the variational method for 
the three-dimensional baffled chamber shown in Fig. 2 is 


2 it r 


/ (0) (u). aJ**) 

(p V -P V ') 22 rdrde 

dz 


= 0 at 


z = z. 


(D-l) 


Again it may be noted that this is simply an equation for conservation of the 
energy flux across the interface between the main chamber and the baffle 
compartments for the rigid-wall case. Corresponding to any pressure gradient 
at the interface, the pressures in the different regions are given by the 

integral equations, i.e.. 


p'“ U/ (r,5,z) = - 


Vl - 1 (M) 


(l u) 

G'‘~ / (r,0,zlr , J e'fZ^) ||t — (r 1 , e^z^r ! dr ’dfl’ 

(D-2a) 


2tt r. 


P (r,a,z) 


c— IS "*• 

■ i f 


G^°^(r#8>z|r , > e’,z l ) - (r', a^z^r'dr'dA* 


(D-2b) 


) 10 ) . 4-V. 

where G and G are the Green s functions in the baffle compartment 
and the main chamber, respectively. The pressures obtained from Eq. D-2a and 
D~2b are not equal at the interface unless the pressure gradient at the inter- 
face used in these equations corresponds to the exact solution of the wave 
equation and the appropriate boundary conditions. 


Because the exact solution is unknown at the outset, an initial form for this 
pressure gradient is assumed and an iteration scheme is then used to improve 
the accuracy of the initial form. This iteration scheme proceeds as follows. 


D-2 



First, some form is assumed for the pressure gradient at the interface: 


3p(^)( £ )(r,0,z 1 ) 

9z 


fu'lfP’) ( o ) ( £ ) 

Then the above integral equations are used to obtai n p'^ A £ ; (r,9 ,z) and p v (r,0,z 
directly. Since, for the exact solution, the pressures must match identically 
at the interface (z = z ), a new approximation for the baffle compartment 
pressure profile at the interface is obtained by equating the pressures: 


(hO(£+i) (o)(0, . v 

P (r,0 > z i ) = p (r,e,z i ) 


(D-4) 


A new pressure gradient is then obtained by differentiating the corresponding 
integral equation and the iteration sequence is complete. The iteration may 
be repeated as many times as desired. 


Assumptions and Definitions 


The assumptions and definitions used to obtain the algebraic forms of the 
iteration and characteristic equations are presented here for reference. 

The £ approximation to the pressure profile on either side of the interface 
(at the interface) is assumed to be represented by an eigenfunction expansion, 
i.e. , 


00(0 


(o)(0 


■ I i <-, 

V 1 V 1 


(^,6,^) - J a (°K f ) + (°) ( r e ) 

0 r 


(D-5a) 


( D- 5b ) 
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where the orthonorml eigenfunctions are 


4^ ( r >0) = J a ^l (<*£^ r)cos(o^ o)/[a^ e^]2 

n Q jll r °V n 0 r,n 0' ,n r n Q> n r n 0 J 


4^ 


(o) 


(r,8) - jJ°> (<, (0) r)cos(or (0) e)/[A (0) « (0) ] 5 


0> m 0 


- 6 '~ r e ' m e Vt m e 

the corresponding eigenvalues are 

(r) (V 1)tt 


= it (n -l) 


\n a (q -0 ) 2 v “0 

6 (JL+1 R 


a Q „ sK- 1 ) 

B s m q e 


th 


a = n~ root of dJ ^ (a^ r )/dr = 0 

r ^ n e^ n T r a e> n e 1 


th 


a.. 


(°) , (0) 




3 n root of dJ (<* ‘ r. )/dr = 0 

r <fy»e 1 


and 


A 


(r) 


n. s n 
0 r 


/ (r) 2 


( (r) r \2 , 

n r i) J 

r,n e ,zl r 


J O 1 ), (r) r 2 

V e a r,n Q> n r r ^ 


A 


a Q »si 
0' r 


, (0) *2 

(a e,m e } 

7 (0) 7^2 

icy / I 

r# 0 ,ta r 1 J 


■ (0) (o) - 2 

Jar (or r, ) 

0,nu r >V m r 


(r) = 

gZ 

f 2 v/r 

n. 

c 

n 

1 Vr 

n l 

(0) 

r 2 v 

m, 

1! 

■1 , 



ll = (1 + 


D-k 


(D-6a) 

(D-6b) 

(D-7a) 

(D-?b) 

(D-7c) 

(D-7d) 

(D- 8 c) 

(D- 8 d) 

(D-8a) 

(D-8b) 



The Kronecker delta is defined by 



m = n 
m ^ n 


(B_9) 


Integrals of products of the orthonormal eigenfunctions are then 



4^ (r,0)^^ ,(r,9)rdrd9 = 

n 0> n r n 9 ^ n r 


6 1 6 t 

n 0' n 0 n r> n r 



4- 


(0) 


m 0^ m r 


(r,e)i|/?^ ,(r,9)rdrd9 = 
m 8 ,in r 


6 

m e ,m e 



in 


/On rH-+l jrl 

I, =/ / ^ (r,e)4^ (r,9)rdrd0 

K'”r) Jq J V m r n 9> n r 

\ n 9, n r / K 


(D-lOa) 


(D-lOb) 


(D-ll) 


The Green’s function for each region has been represented by an eigenfunction 
expansion also, 


G^(r|r') s - V V 4 /^ (r,©)^^ (r‘, e’)F^ 

n^l nt=L V n r V n r V n i 


(zlz') (D-12a) 


9 r 


G (0) (?|?') = + 


Z 2^ 0) m ^ r, ^ 0, ) F i o) m (z i z!) (D - i2b 


m =1 m =1 m 0 ;D V ’ m 0 ,m r 

U T 


m Q .m 
Q J r 
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where 


f ( h) ( z j z * ) 2 < 

n_,n 


cos 


(z-z ) COS [^2^2 n ( z p“ z *) + n ] 

' Q’ r L z > n o> n vd 2 Rz ^ n A^ n r J 


B’ r 


e ; r ' 


~m z&r~ 

a n 

z 'V n r n 0; n r 


z < 


cos [ar^ (z-z )]cosfo'^^ (z ? -z) + ] 

L z^n/ l'J L z^n/ 2 Rz » n Q» n r J 


(7) T(mT“ 

a £2 

z,n Q ,n r n 0 ,n r 


z > 


V m r 


(z|z*) s i 


cos 


[ar(°) (z-zJ + &(°^ IcosT^ 0 ) ( z l-z')l 

L ,lO ° J K rr, rr, -1 L J 


z ,V m r 


( 0 ) ( 0 ) 

a £2 

z^mQ^nij, m 


cos [ar^ (z'-z ) + 

L z,m 0 ,m r v o' L 


’ z > m Qj m r 


cos 


[*[°l ^ z i" z )] 

1 z ; m e ^ r 1 J 


'^r— (or - 

z,m 0 ,in r m 0 ,m r 


in which 

(^) 


c = - ( a i^ " f 

z } n e ,n r 


r > n d> n r 


a 


( 0 ) 


z 


5 Jk 2 - (a (0) )2 


£2 


(M-) 

n 0> n i 


- sin 


in[ar^ Ug-z ) + ] 

z^n e ,n r * i z,n flJ n r J 


fl^ u r 


£2 


( 0 ) 


m 0> m r 


= sin R°m ™ (z r z o ) + i (0) 


z > m 0^ m r 


and the finite admittance boundary conditions give 


- (m ) 

I 

x z,HQ,n r 


-1, 


» 


“ tan (ikp 7 j or * ) 

z 2- z,n 0 ,n r ; 


»(o) 

"*2* ^ 133.^ ^ XIL^, 


- tan" 1 (ik0 _ ) 

z o+ ^“^“r 


(D-13a) 


(D-13b) 
> z ! 

(D-l4a) 
(D-l4b ) 
(D-l4c) 
(D-l4d) 

(D-15a) 

(D-15b) 
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Finally, 


%; n r 


(z) 


( , sinfc/^ (z -z) + 

= + >) L£'_^P_r... 2 Rz ' n 0' n r 

z ' n 0> n r r (l-0 ( v (n) 

cos a (z 0 -z, ) + £ 

L z,n ,n 2 1 *z,n ,a 


d’ r 


(D-l6a) 


i l(0) (o) sin far' (z-z ) + & ^ j 

H ' (z) 3 - a { } L z,m 0 ,m ^ o' bz,m 0 ,m J 

^ r 0 r r (0) (0) 1 

a V ' (z-z) + ^ ‘ ! 

L z,m 0 ,m T . v 1 o' ^z,m 0 ,m r J 


cos 


(D-l6b) 


B 


(0)02) 


V m r 


■ms cxi-v- ■ 


i.°) 


M=1 i n 0 =l n r =l 


m 0 ,mr 

V n r 
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Final Equations 


When series expressions are used to represent the interface pressure profiles 
in both the main chamber and the baffle compartments, the iteration scheme 
proceeds as follows : 


MU) 


I X ^ ) i* ) * B w n ('.») 

J. £—4 r Z—s HqjTi j, 


n 0 =l n r =l 


(D-l8a) 


3P 


(M)(«) 


3z 


8 s y < } ‘ e) d t^* 00 mo) 

x ^ n 0 > n r n eJ n r x n„,ru. 


n =1 iy=l 


(B-lSb ) 


r ‘■ L Q» u -r 


( 0 )( 2 ), . _ 

P (r, 0 ,z, ) = 


B 


( 0 )(«) 


r S I u( ' o) , . 

®fl = l m =1 Z 1 


m 0 ,,m r , (0 ) , < _ v" ir-~ ■ (0)(£) (o) 

r K „ (r,0) = \ y a f ; (r,0) 


"0 _x V x 


‘0 ,m r 


I I 

m 0 =l m =1 
0 r 


^“r 

(B~l8c) 
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dP ( 


\v s Q } z 


>. y y B (o)(s v (o) (r,.j. y y (a (o)w w< 0) <*,» 

m^l m -1 V m r ”«’”V „4, 

0 r 


m =1 m =1 
q r 


^ (0) (r, 0 ) 

m„,m 7 
0 7 r 


0 7 r 
(D-l 8 d) 


(m)(C+ 1) 


‘-^2 * \\{\ I . 

n^=l n r =l 17 x 17 x n 0 =1 n r 1 \ m e =1 “V 1 - 


V (M) (M> 

V m r (V m r)jv n r 

n Q ,ri T 

With the same series expression representations, the general characteristic 
equation reduces to an algebraic equation (which can be solved numerically) 


(D-l8e) 




m 0 =l m =1 

u x 


6 r 


00(«) 


n 0 7n r 


V M) (o' 

Rn „’ n r 1 ) 


= 0 
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COMPUTER PROGRAMS 


The computer programs used to solve for the fundamental acoustic modes of a 
three-dimensional baffled chamber with gain/loss end-wall boundary conditions 
(see Figure 2) are described in this section. All programs are written in 
GE-400 Series Fortran* using complex arithmetic. The function of the com- 
ponent programs, the programs themselves, and the relationship between the 
functions, the equations and the lines in the program are detailed below. 

The computer logic is illustrated by the flow chart shown in Fig. 14- 

*C?B-l473, GE-toO Series Time-Sharing FORTRAN manual. General Electric Company, 
Information Systems Division, 13^30 North Black Canyon Highway, Phoenix, 

Arizona, 1970. 
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RN5353 - INPUT FILE 


This file lists data for input both to the r integration program and to 


main program. 




Specification 

of Numbers of 

Terms : Line 

0100 

MTOVMX = 

(m ) =5 

0'max 



MROVMX s 

(m ) =3 

r max 



NTUVMX = 

^ n 0^max - ^ 



HRUVMX s 

^ n r "max ~ ^ 



Specification 

of Integration Parameters : 

Line 

MINGRL = 

80 



EPSSER = 

10~ 25 



IPRDET = 

0 



Specification 

of Zeros of J 1 

! and Values 

of J: 


ATOM(l) = J e °l = 0.00 

AROM(l.l) = ( J° ) r, ) = 0.00000 . J (o)i a ^°} ,r. ) = +1.00000 
r.1,1 1 5 a K J r,l,l 1' 

0,1 

AR0M(1,2) s (a^ 0) r ) = 3-83171, J (o)^^ p r ) = -0.40276 

Ty±-)C. X i^X^X x 


AT0M( 


&EOM(l,3) 2 (7°? 3 rj_) - 7-01559, j (0)(“r°! jlp = + 0 . 30012 

2 ) , = j -’ 


6,2 


.00 


AR0M(2,1) = (a^ r.) = 1.84118. J (o )(^°' > r ) = +O.58187 

r,2.1 1 3 cr ' r,2,l 1 

J ' 9,2 
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RI5353 - INPUT GENERATION PROGRAM -r INTEGRATION 

This program is used to calculate the integrals of products of Bessel func- 
tions of different orders and to create a file of integrals for use by 
the main program. 


a. Input: Lines 1120, 1150 , 1170, 1190: From File RN5353 

MTOVMX = ( m g) max - if terms, main chamber e series expansion 
MROVMX = (n^ ) max = =§= terms, main chamber r series expansion 
NTUVMX = (n ) = # terms, baffle compartment 6 series expansion 

NRUVMX s ( n r ) max = H "terms, baffle compartment r series expansion 


MINGRL = ( m ) max = if terms, r integral series approximation 

EPSSER e e = term value, r integral series approximation cutoff 

IPRDET a Inactive 


ATOM(MTOD) 

arom(mtod,mrod) 

bessom(mtod,mrod) 



atum(ntud) = 
arum(ntud,nrud) s 
bessum(ntud,nrud) 3 


a 


a 


,(l0 

e,n 


J (n) (a^ r 1 ) 


e r n e 


D-ll 



Perform r Integration: Lines 1210- 1780 + Subroutines 


FINHM(MTOD,MROD,NTUD,NRUD) = £ 


(g) 


fVM 
'VV 


(Eqn. D-ll) 


SUBR0UTINE SERIEF - Calculates series approximation to the 

* 

r integral 




J M (ar)J v (br)rdr = 


_ M , v 
2/arf ,br, 

r 'T' 


r(v+i) 


^ (-l) m (^) 2m p F 1 (-m,-M-m;v+l;b 2 /a 2 )] 


2 1 

( 2+M+V+ 2m )m.' r(jn + m+ 1 ) 


where 

T is the Gamma Function - SUBR0UTINE GAMMAF 

and 


p F_ (or,0;p;z) 

c - 1 


91 (a) (Q) 

V „ n 

Z n.' ( 


n=0 


(P) 


n 


(«) 

(«), 


n 


a(a+l)(ar+2). . . (a+n-l) 
1 


is a Hypergeometric function’ 


that is evaluated by SUBROUTINE HYPERF 


o ® 


Output: Lines 1790.? 1810 

FIMEM filed in RT5353 

MTSRMX = ( m ) max = if terms actually used in series approx. 

EPSSER = e 


printed 


"Watson, G. N« , A Treatise on the Theory of Besael Functions , Cambridge 
University Press, London, 1962, pg. 148. 

4, 

! Ibid, p. 100 
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TIT355 - INPUT GENERATION PROGRAM 


This program is used to generate a file of theta integrals needed by the 
main program. 


a. Input: Line 10 4-0 

THTO = e Q 
MUBAR = M 

MTOVMX = (m ) = # terms, main chamber 0 series expansion 

® max 

MTUVMX = ( n Q) max = $ terms, baffle compartment 0 series expansion 


b. Perform 0 Integration: Lines 1050-14-70: (Eqn. D-ll) 


FINTTM( MUD, MTOD, NTUD ) 






c. Output: Lines - 1480- 1510 


FINTTM filed in T0T355 


THTO 

MUBAR 

MTOVMX 

NTUVMX 


printed 
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B3D353 - MAIN PROGRAM 


This program is used to obtain the eigenvalue (frequency and damping coefficient) 
as veil as the series expansion coefficients for the interface pressures and 
axial velocities for a particular three-dimensional baffled chamber with 
gain/loss end-wall boundary conditions. Complex quantities are noted. A full 
sample calculation is given at the end of this Appendix. 


Input Files Needed 

RN5353 - Bessel Function Data 


RT5353 - r Integrals £ 


Integrals g 


(°h 

(Vt 

V n r 

(2) 

( m e) 


(Eqn. D-ll) 


(Eqn. D-ll) 


Input: Lines 1290 , 1300 - 1310 , 1320 - 1330 , 1340 , 1350 

= ( m e^ m ax s ^ terms, main chamber 0 series expansion 
= ( m r ) max = # terms, main chamber r series expansion 
NTUVMX 5 s #terms, baffle compartment 9 series expansion 

U IBclX 


NRUVMX 


In 

Input 

Data 

Files 


( n r ) m ax 5 ^ terms, baffle compartment r series expansion ' 


IDJM1 

DUM1 

IDJM2 


Inactive (dummy variables) 


ATOM(MTOD) 

arom(mtod,mrod) 


- (0) 

a e,m e 

s r-L 


BESS0M(MT0D,MR0D) = J (0)(a^ 0 ^ 


r ) 


r,m 0 ,m r ! 


ATUM(MTUD) 

aeum(mtud,mrud) 


si a 




5 <4,n r r l 


BESSUM(MTUD,MRUD) s J (m-) ( a (^) r, ) 

an - r.n " 


0,n, 


0 ,n r 
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b. Input 1: Line 1370 


MUBAR 


R1 

Z0 

Z1 

Z2 

BETAOP 


= z 


= z_ 


= z. 


= p. 


"o+ 


BETA2L = S 


(Complex) 

(Complex) 


Definition: Lines 1400-1420 


BSUM01 = = +0.5 


- 1.0 


BSUM02 = b^ 2 ^°) 

0,0 

BSUM03 = b^ 2 ^°) = +0.5 

0,0 


initial pressure distribution 
(other coefficients are zero) 


d. Input 3: Line 1440 


MTOVMC = # 6 terms retained, main chamber series 

MROVMC = # r terms retained, main chamber series 

NTUVMC = § d terms retained, baffle compartment series 

NRUVMC = # r terms retained, baffle compartment series 

JIT = number of iterations performed 


e. Input 4: Line 1460 


ITYPE = 


0 for step calculation of characteristic function 

1 for root finding calculation of characteristic function 
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f. Input 5: Line 1500 (ITYPE = 0) 


EIR1ST = kr^ initial guess (Complex) 

EIR1SP = kr^ increment (Complex) 

NEIR1S = number of calculations 

| 0 for simple printout of characteristic function 


XDETAL = 


1 for detailed printout of series expansion coefficients, 
pressure and velocity profiles 


g* 


Input: Line 1530: 

IPBSUM ) 


( ITYPE = 0 ) ( IDETAL = l) 

f BSUM(MUD,NTUD,NRUD) 


IPVBUM I 


0 for no print 

1 for printout 


of< 


VBUM( MUD, NTUD, NRUD) 


IPASOM 

IPVAOM 


AS0M(MT0D,MR0D) 
k VA0M(MT0D,MR0D) 


b M(i) 

V n r 

b (^)(i)|_)(^) ( j 

n 0 > “r n e 3 n r 

a (o)C4) 

V m r 


m 0> m r 1 


d* r 


Input: Line 1550: 


MPBSUM 

MPVBUM 

MPASOM 

MPVAOM 






0 

1 


(ITYPE = 0 ) ( IDETAL = l) 


for no print 
for printout 


f 

P 


3P 


of < 


P 


dP 


(f*)(i) 

(K)(i) 

(o)U) 

( o )(0 


(r^z.^ 
(r,0,z 1 )/az 
(r, 0,z^) 

(r, 6,z^)/dz 


Input 6: Line I 58 O: (ITYPE = l) 

EIB1ST = kr^ initial guess (Complex) 

EIB1FD = kr^ finite difference (Complex) 

ME HUM = number of calculations 

CHFNTL = tolerance -final characteristic function value 
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i. Definition: Lines 1700-2110 (Eqns. D-7a,b,c,d;D-l4a,b;D-15a,b;D-l6a,b) 


j • 


JI 

= 0 

(Complex) 

AZOV 

z,m 0 ,m r 

( Complex ) 

DELLOV 

= |(°) 

z,mg,mr 

( Complex ) 

HLOM ( MTOD , MROD ) 

= tf 0) (O 

(Complex) 

AZUV 

z,ng,n r 

(Complex) 

DELLUV 

z ’V n r 

( Complex ) 

HRUM( NTUD, NRUD) 

= g M U ) 
N^V l' 

( Complex ) 

BSUM( MUD , NTUD , NRUD ) 

- b OO(0) 

n 6’ n r 

(Complex) 


Calculation - Iteration Dependent Variables: Lines 2140-2380 (Eqns. D-17, D-l8) 


VBUM(MUD,NTUD,NRUD) = b (fX)(C) td (,X) (z ) (Complex) 

n 6 }XX T V n r 1 


BL0M( MTOD, MROD) 

= B (o)(s) 

m 0 ,mr 

( Complex ) 

vaom(mtod,mrod) 

. (0 ) ( C ) 

- a 

p in 

= a (0)(«)y( 0 ) ( 

m e }Ta r m e> m T 1 

(Complex) 

asom(mtod,mrod) 

( Complex ) 


k. Perform Iteration: Lines 2390-2590: (JI < JIT): (Eqn. D-l8e) 


JI = JI + 1 
BSUM( MUD, NTUD, NRUD ) 


(H)(£+l) 

b 

n 0' n r 


(Complex) 
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Compute Characteristic Function: Lines 2620-2850: (Eqn. D- 19 ) 



(Complex) 


(Complex) 


CHAEFN - SUM1 - SUM2 (Complex) 
General Printout: Line 2890 

1IGMR1 = kr x (Complex) 
SUM1 (Complex) 
CHAEFN (Complex) 
SUM2 ( Complex ) 


'Detailed Printout: Lines 2920-3140: (IDETAL = l) 

(According to inputs g, h) 

Calls SUBROUTINES f PV 3 DBC \ for interface pressure/axial- velocity profiles 

TV3DMC ; 

Increment Eigenvalue: Lines 31^0-3^80 

ITYPE = 0 EIGNR1 = EIGNR1 + EIR1SP (Complex) 

ITYPE = 1 =2> EIGNR1 = EIGNR1 - CHARFN 

d( CHARFN) 
d(EIGNRl) 


Return Options: 

Lines 3500-3540 

IG0T0 = 1 

Input 1 

1G0I0 = 3 

Input 3 

IG0T0 » 4 

Input 4 

IGjZSl^ = 5 

Input 5 

IG0T0 = 6 

Input 6 



SUBROUTINES CALLED BY MAIN PROGRAM 


PV3DBC - Baffle Compartment Interface Pres sure/ Axial Velocity Profiles 

This program takes coefficients of pres sure /axial- velocity series expres- 
sions and. generates the appropriate radial and tangential profiles at 
the interface on the baffle compartment side. 

a. Input: Line 5300 

NRPNT s number of equally spaced r printout points 
TVAL s Q value for r profile 

NTPNT = number of equally spaced 6 printout points 

RVAL s r value for 6 profile 

TST0P = last 6 value printed 

(NRPNT s NTPNT s 0 for return to main program) 

b. Output: Line 5520: (Eqn. D-6a) 

Value of pressure/axial velocity at r location 
r location 

c. Output: Line 5970: (Eqn. D-6a) 

Value of pressure/axial velocity at 6 location 
9 location 
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PV3DMC ~ Main Chamber Interface Pres sure/ Axial Velocity Profiles 

This program takes coefficients of pres sure/ axial velocity series expres- 
sions and generates the appropriate radial and tangential profiles at the 
interface on the main chamber side. 

a. Input: Line 6280 

NBPNT = number of equally spaced r printout points 
TVAL = 6 value for r profile 

NTPNT = number of equally spaced 9 printout points 

RVAL = r value for 6 profile 

TST0P = last 6 value printed 

(NRPNT = NTPNT = 0 for return to main program) 

b. Output: Line 6530: (Eqn. D-6b) 

Value of pressure/axial velocity at r location 

r location 

c. Output: Line 6860: (Eqn. D-6b) 

Value of pres sure/ axial velocity at 6 location 

9 location 
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AUXILIARY SUBROUTINES 


BESFNP - Bessel Function Calculation Program 

This program takes the order and argument from the variable list and, for 
Bessel functions of integer or half-integer order, returns the value of 
the Bessel function using standard recursion formulae. 

a. Input: Line J000 (Variable List) 

ORDER = order of Bessel function (integer or half integer) 

ARGMNT = argument of Bessel function 

b. Output: Line 7000 (Variable List) 

BESFNV s value of Bessel function J 0 Rr)ER (ARGMNT) 

BESSOP - Zero Order Bessel Function Evaluation 

This program takes the argument from the variable list, upon call by 
BESFNP, and returns the value of the zero order Bessel function of that 
argument. 

a. Input: Line 8000 (Variable List) 

ARGMNT = argument of Bessel function 

b. Output: Line 8000 (Variable List) 

BESSOV 5 value of zero order Bessel function J (ARGMNT) 
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BESSIP - Order One Bessel Function Evaluation 


This program takes the argument from the variable list, upon call by 
BESFNP, and returns the value of the order one Bessel function of that 

argument. 

a. Input: Line 9000 (Variable Test) 

ARGMNT s argument of Bessel function 

b. Output: Line 9000 (Variable List) 

BISSIV s value of order one Bessel function J^( ARGMNT) 
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1 


RN 5353 



0100 

3 

3 

5 3 

0200 

80 

1 . 0 E -25 0 

1010 
10 1 1 

0-00 

0.00000 

+ 1 .00000 

10 1 2 


3.831 71 

- 0 . 40276 

101 3 


7.01 539 

+ 0.30012 

1020 

1021 

1 .00 

1.84118 

+ 0.58187 

1022 


5.33144 

- 0-34613 

1023 


8 .53632 

+ 0.27330 

1030 

1031 

2.00 

3-05424 

+ 0 • 48650 

1032 


6.70613 

- 0-31353 

1033 


9.96947 

+ 0-25474 

10 40 
1041 

3-00 

4*20119 

+ 0 . 43439 

1042 


8 .0 1 324 

- 0.29116 

1043 


1 1 -34592 

+ 0-24074 

1050 
10 31 

4.00 

5.31755 

+ 0.39965 

1032 


9 .28240 

- 0.27438 

1033 


12.68191 

+ 0.22959 

2010 
201 1 

0-00 

0 .00000 

+ 1 .00000 

2012 


3.83171 

- 0 • 40276 

201 3 


7.01 559 

+ 0.30012 

2020 

2021 

1 • 30 

2.46054 

+ 0 • 52534 

2022 


6 .02929 

- 0.32806 

2023 


9 . 261 40 

+ 0.26330 

20 30 
2031 

3.00 

4.201 19 

+ 0 . 43439 

20 32 


8.01 524 

- 0.29116 

2033 


1 1 .34592 

+ 0.24074 

20 40 
20 41 

4.50 

5.86842 

+ 0.38601 

20 42 


9 .90431 

- 0.26739 

20 43 


1 3.33793 

+ 0.22479 

20 30 
20 31 

6.00 

7.50127 

+ 0.3541 4 

20 32 


1 1 . 73494 

- 0.2501 7 

2033 


1 5.26818 

+ 0.21261 
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1 


R I 5 3 5 3 


1000 DIMENSION ATOM C 5 ) * AROM (5j3)> BESSOiM C 5* 3 ) 

10 1 0& * ATOM (5)* ARUM C 5* 3 ) * BESSUM C b* 3 ) 

1Q20& * F I NTEM C 5* 3* 5* 3 ) 

1030 COMMON A TOO* AROV*ATUV* ARUV* AROV02* ARUV02* A TO VP bATUVPli RARAR2 
1040&*MINGRL* EP3SER* I PRDET* MTERM* RM* SER I EV* HYPERV 
1050 DATA MTOVST* MROVST* NTUVST* NRUVST/ 1> N 1 j 1/ 

1060 CALL 0PENF ( 1 * "RN5353") 

1070 CALL 0PENFC2* "RT5353") 

1080 2 F 0RMAT (415) 

1090 A FORMAT ( I 5, El 0 .2* I 5) 

1100 6 F0RMATCF5.2* 3C/5X* 2F9.5) ) 

1110 12 F0RMATC 1X 5 1P5E1 2.6) 

1120 READC 1 * 2) MTOVMX* MROVMX* NTUVMX* NRUVMX 
1 1 30 P R I N T $ *'MTOVMX="*MTOVMX* " MROVMX="*MRO VMX 
1140 PR I NT* "NTUVMX 2 "* NTUVMX » " NRUVMX="* NRUVMX 
1150 READC 1*4) MINGRL* EPSSEft* IPRDET 
1160 PRINT* ”MINGRL="*MINGRL*” EPSSER="* EPSSER 

1 1 7 0 READC 1*6) ( ATOM C MT ) * C AROM CMT* MR) * BESSOM C MT* MR) * MR= 1 *MROVMX)* 

1 1 SO&M T= 1 * M TO VMX ) 

1190 READC 1*6) C ATUMC NT) * C ARUM C NT* NR) * BESSUM C NT* NR ) * NR- 1 * NRUVMX) * 

12007 N T = 1 * NTUVMX ) 

1210 M T 0 V N D = M T 0 V M X 
1 2 2 0 Mi h 0 VND=MR0VMX 
1230 N T U V N D = N T U V M X 
1240 N R U V N D = N R U V M X 
12^0 MT£RMX=0 

1260 DO 90 MTOV=MTOVST*MTOVND 
1270 D0 80 MROV=MROVST*MROVND 
1280 D0 70 NTUV=NTUVST* NTUVND 
1290 D0 60 NRUV=NRUVST* NRUVND 

1300 ATQV=AT0MCMT0V) 

1310 A TO VP 1 =ATO V+ 1 .0 

1320 A R 0 9 = A R 0 M C M T 0 V * M R 0 V ) 

1330 AROV02=AROV/2.O 
1340 ATUV= ATUM ( NTUV ) 

1350 ATUVP1 =ATUV+1 .0 
1360 ARUV=ARUMCNTUV*NRUV) 

1370 ARUV02=ARUV/2.O 

1380 IF C CAROV.EQ .0 .0) .0R. CARUV.EQ.O .0) ) G0 T0 10 

1390 AROATO=AROV02**ATOV 
1 40 0 RA TARO = < ATO V /A RO V ) **2 
1410 ARUATU=ARUV02**ATUV 
1420 RATARU= C ATUV/ARUV ) **2 
1430 RARAR2=CARUV/AR0V)**2 
1440 GO T0 40 

1450 10 IF C CAROV «EQ .0 .0) .AND. CARUV.EQ .0 .0) ) G0 T0 30 

1460 IF CARUV.EQ.O.O) G0 T0 20 

1470 AIQ V=ATUV 
1480 ATOVP 1 =ATUVP 1 
1490 ATUV=Q • 0 
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2 


RI5353 CONTINUED 


1500 ATUVP 1 = 1 .0 

1510 AROV=ARUV 

1520 AftOV02=AftUV02 

1530 ARUV=0 . 0 

1540 ARUV02=0 • 0 

1550 20 AR0AT0=AR0V02**AT0V 

1560 RATAR0=CAT0V/AR0V)**2 

1570 ARUATU= 1 .0 

1580 RATARU=0 • 0 

1590 RARAR2=0 • 0 

1600 G0 TO 40 

1610 30 AR0AT0= 1 .0 

1620 RATAR0=0.0 

1630 ARUATU= 1 .0 

1640 RATARU = 0 • 0 

1650 RARAR2= 1 .0 

1660 40 CONTINUE 

1670 CALL SERI EF 

1680 CALL GAMMAF C ATUVP 1 * ATUVFT) 

1690 FINTEV=(AROATO*ARUATU*SERI EV) /ATUVFT 
1700 BESS0V = BESS0M (MTOV* MROV ) 

1710 BESSUV=BESSUM ( NTUV* NRUV) 

1720 REF ERV = 0 . 5*SQRT< < 1 -O-RATAROl+C 1 .O-RATARU) ) *AB S C BESSOV*BES SUV) 

1730 FI NTEM ( MTOV* MROV* NTUV* NRUV ) =F I NTEV/REFERV 

1740 IFCNTERM.GT.MTERMX) MTERMX=MTERM 

1750 60 CONTINUE 

1760 70 CONTINUE 

1770 80 CONTINUE 

1780 90 CONTINUE 

1790 WRI TEC 2* 12) FINTEM 

1800 CALL CL0SEFC2) 

1810 PRINT* t 2* m MINGRL = "*MTERMX.» " EPSSER=’’* EP SSER 
1820 STOP 
1830 END 
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2000 SUBROUTINE SERI EF 

20 1 0 COMMON ATOV* AROV> ATUVj ARUV* AROV02.» ARUV02* ATOVP 1 * ATUVP I * RARAR2 
20 20 & * M I N G RL > EPSSER* IPRDET-. MTERM, RM, SERI EV> HYPERV 

2030 ARO022=ARO V02**2 
20 40 ATPAT=ATOV+ATUV 

20 50 CALL GAMMAF ( ATOVP 1 .» ATOVFT) 

20 60 MTER21 = 0 

2070 TM= 1 .0/C(ATPAT+2.0)*ATOVFT) 

2080 CALL HYPERF 
2090 SERI £V=TM*HYPERV 

2 1 U 0 H i v i = 0 • 0 

2110 DO 10 MTERM= 1 >MINGRL 

2120 RiM = RM+ 1 • 0 

2130 RMT2P=2 •0*RM+ATPAT 

21 40 TM=- CARO 022/ CRM* (RM+ATOV) ) > * C RMT2P/ ( RMT2P+2 . 0 ) >*TM 

2150 CALL HYPERF 

21 60 SERIEV=SERIEV+TM*HYPERV 

21 70 IF CABSCTM+HYPERV) .LT.EPSSER) RETURN 

2180 10 CONTINUE 
2190 RETURN 
2200 END 



1 


HYPERF 


3000 SUBROUTINE HYPERF 

3010 COMMON ATOVjAROVjATUV^ARUVjAROV02jARUV02jATOVP 1 , ATUVP1 .» RARAR2 

3020& .* M I NGRLj EPSSER* IPRDET* MTERM, RM, SERI EV, HYP ERV 

3030 RMP 1 =RM+ 1 • 0 

30 40 TN= 1 .0 

30 bO HYPERV=TN 

3060 IF ( MTERM . EQ • 0 ) RETURN 

3070 RN=0 • 0 

3080 DO 10 N= 1 , MTERM 

3090 RN= RN+ 1 • 0 

3100 RMP 1 LR = RMP 1 - RN 

3110 TN= C ( RMP 1 LR* ( RMP 1 LR+ATOV ) ) /C RN*C RN+ATUV) ) ) *C RARAR2*TN) 

3120 10 HYPERV=HYPERV+TN 
3130 RETURN 
3140 END 
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GAMMA F 


40 0 0 SUBROUTINE GAMMAF C VI N, V0UT ) 

4010 DV I N = V I N~ 1 .0 
40 20 I N = DV I N 
40 30 D V I N T = I N 

40 40 IF ( ( DVI N-DVINT) • GT.O . 1 ) G0 T0 20 

40 SO V 0 U T = 1 * 0 

4060 10 IF CDVIN.LT. 1 • 1 ) RETURN 

4070 V0UT=DVIN*V0UT 
40 8 0 DV I N = DV I N- 1 .0 

40 9 0 G0 T0 10 

4100 20 V0UT= 1 • 7724339 

4110 25 IF CDVIN.LT. 0.4) RETURN 

4120 V0UT=DVIN*V0UT 

4130 DV I N = DVI N- 1 • 0 

4140 G0 T 0 25 

41 50 END 
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READY 

EDIT WEAVE RI 53 53, S ERI EF* HYPERF, GAMMA F 

READY 
0001 SN DM 

RUN 


generate 3TS3$3 


RI 5353 


20:03 NAR JAN. 12,1971 


MTO VMX= 

5 MRO VMX= 

3 


NTUVMX= 

5 N R U V M X = 

3 


MI NGRL= 

80 EPSSER= 

1 • 000000000 E- 

25 

MI NGRL= 
STOP 

47 EPSSER= 

1 . 000000000 E- 

25 

RUNNING TIME: 

56.2 SECS 

1/0 TIME : 

1 • 


2 SECS 
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1 


RT53S3 


0000 1 
00002 
0000 3 
0000 4 
0000 b 
00006 - 

0000 7 
00008 
GO 00 9 
000 10 

0001 i 
000 1 2 
0001 3 

000 1 4 

0001 b 
00016 - 
000 1 7 
00018 •“ 
000 1 9 
00020 
00021 
00022 

00023 

00024 - 

00025 

00026 

00027 - 

00028 

00029 

00030 - 

00031 

00032 - 

00033 
0003 4 
00035 
000 36 
0003 7 

00038 

00039 
000 40 
000 41 
00042 * 
000 43 

00044 

00045 - 


1 •QOOOOE+OO 
2 . 10339E-06 
9 . 49233E-07 
9 . 44745E-01 
3 .27231 E-01 - 
1*71 789 E- 02- 
8 » 8 35 53E~0 1 
4 • 620 54E-0 1 - 

7 • 6 4 6 7 2 E - 0 2 

8 3 4738 E- 0 1 
5.24798E-01 - 
1 .6 48 50 E-01 

7 . 9 S438E-0 1 
5 » 55671 E-01 - 
2-341 42E-01 
2. 1 0339E-06- 
9 .99997E-01 
4* 6 1 341 E-0 7 
2.66278E-01 
7.98422E-01 
5.371 75E-01 - 
3.5661 IE-01 
5.60109E-01 
7 • 32880E-0 1 - 
3*97926 E-01 
3 . 77426E-01 
7.62856E-01- 
4.19145E-01 
2-41 269 E-0 1 
7.24129E-01- 

9 * 49233E-07- 
4.70261 E-07- 
9 .99972E-01 

1 » 3341 6 E-0 1 

3 » 6 5832E-0 1 
6 • 72104E-01 
1 .98336E-01 
4.27188E-01 
2.61464E-01 
2 .35841 E-01 

4 03 5 60 E-01 
3 « 30600E-02 
2.59524E-01 
3*55163 E-01 
2.22725E-01 


9 .6731 3E-01 
2*1 1320E-01 
1 .0041 5E-01 
9.96511 E-01 
8 .2021 5E-02 
1 * 30458E-02 
9 .67401 E-01 
2.52876E-01' 
1 .3 48 42E-02 
9 . 32908 E-0 1 
3 • 52624E-0 1 
7.221 4 3 E - 0 2 
9 .00 522E-0 1 
4. 1 3531 E-0 1 
1 .31 372E-01 
2.501 42E-01 
8 • 8 4508 E-0 1 
2 • 99353E-0 1 
7.89033E-02 
9 .851 39E-01 
1 • 4729 1 E-0 1 
2.29840E-01 
8 . 54296E-01 
4.66097E-01 
3.1 4024E-01 
6.99101 E-01 
6 • 25355E-0 1 
3 * 65685E-0 1 
5 • 59 1 26E-0 1 
6*86021 E-01 
3.73455E-02 
4.04626E-01 

8 * 1 7943E-01 
2.27360E-02 
1 * 36086E-0 1 

9 * 72363E-0 1 
8*61 584E-02 
3*7 4899 E-0 1 
7.1 7394E-01 
1 * 36 730 E-0 1 
4*671 58E-01 
4. 16991 E-01 
1 .75743E-01 
4.S9015E-01 
1 » 62892E-0 1 


9 • 229 28 E-0 1 
3 • 0 5500 E-01 
1 • 596 1 2E-0 1 
9 • 9 76 57E-0 1 
6 • 593 54E-02 
1 .56329E-02 
9.941 50 E-0 1 
1 • 0 79 63E-0 1 - 
3.72199E-03- 
9 . 74331 E-01 
2.23533E-01 - 
2. 57870E-02 
9 . 50839 E~01 
3.01 520E-01 - 
6.93630E-02 
3.84780E-01 - 
7.1 3454E-01 
4.03202E-01 
6.76728E-02- 
9 .89372E-01 
1 . 19312E-01 
1 .04664E-01 - 
9 . 71 726E-01 
2. 1031 4E-01 
2.12267E-01 
8 . 74342E-01 
4.34598E-01 - 
2.84323E-01 
7 . 60339 E-01 
5.6 5383E-0 1 - 
1 • 1 7572E-02 
6 • 30548E-0 1 - 
5.25268E-01 
8 .36204E-03 
1 .26406E-01- 
9.7901 7E-01 
2 . 37390E-02- 
1 . 9 4447E-0 1 
9 • 40930E-0 1 
6 . 40396E-02 
3 • 748 79 E-0 1 
7 . 34236E-0 1 
1 .02068E-01 
4.681 26E-01 
4*9 7724E-0 1 


8 -83553E-0 1 
3.5661 1 E-01 
1 .98336E-0 1 
9 .84623E-0 1 
1 • 6 49 9 9 E- 0 1 
4« 90277 E -02 
1 .00002E+00 
1 • 40526E-06 
1 . 4746SE-06 
9 .92687E-01 
1 .20 5 78E-01- 
5.01 536E-03 
9 .7761 7E-01 
2.0 78 1 7E- Or 
3*2 46 78 E-0 2 
4.62054E-01' 
5 » 60 1 09 E- 0 1 
4.271 88E-01 
1 . 74684E-01 
9 .28085E-01 
2 • 90 780 E-0 1 
1 * 40 736E- 06 

1 .00000 E+00 
1 . 46999E-05 
1.181 18E-01 
9 . 623 48 E-0 1 

2 * 44656E-0 1 
2 .022 1 0 E-0 1 
8 .84344E-01 
4.1671 IE-01 
7.64672E-02 
7 * 32880E-0 1 
2.61 456E-01 
1 • 1 9396E-03 
3.33400E-01 
8 • 5441 6E-01 
1 .25647E-06 
6 . 78837E-06 
1 .00003E+00 
2.22356E-02 
2 • 28 593E-0 1 
9 . 16220E-01 
5.1 23 79 E- 02 
3 • 74572E-0 1 
7.4371 3E- 0 1 


8 . 49802E-0 1 

3.8 720 6 E-0 1 
2 . 253 1 2 E-0 1 
9 • 6 7039 E- 0 1 
2. 3 52 10 E-01 
8 • 1 1 803E-02 
9 .96421 E-0 1 
8 . 33344E-02 
1 .3571 7 E - 0 2 

9.9 9 329 E-0 1 
3. 6 46 10 E-0 2 
7.47042E-04 
9 .9 1 701 E-01 
1 • 280 1 3E-0 1 
1 . 191 65E-02 
5.08659E-01 
4- 32444E-0 1 
4. 1 58 70 E-01 
2 • 53002E-Q 1 
8 « 44564E-0 1 
3.96792E-01 
8.46711E-02 
9.81 608 E-0 1 
1 . 62905E-0 1 
3.62867E-02 
9 .96461 E-01 
7.51 745E-02 

1 . 26 580E-0 1 
9 . 55530E-0 1 

2 • 6 59 40 E-0 1 
1 .37 432E-0 1 
7 • 63758E-0 1 
5. 1 9486E-02 
2.84210E-02 
4. 78956E-01 
6 . 83296E-0 1 
2.65251 E-04 
1 .71 199E-01 
9 . 588 58 E-01 
3.31 560E-03 
7.371 66E-02 
9 .898 49 E-0 1 
2.04375E-02 
2.51720E-01 
8.9381 5E-01 
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TI T3 d b 


1000 DIMENSION F I NTTM C 3* 5, 5 ) 

1010 CALL BPENF C 1 , "T0T355") 

1020 2 F0RMATC IX, 1P3E1 5*8) 

1030 PRINT," INPUT: THTO , MUBAR, MTOVMX, NTUVMX" 

1040 INPUT, THTO, MUBAR, MTOVMX, NTUVMX 
1050 PI=3. 1 41 5926536 
1060 FMU8AR= FLOAT (MUBAR) 

1070 THTDIF=2.0*PI /FMUBAR 
1080 PTHTDF=PI*THTDIF 
1090 00 1000 MUD=1, MUBAR 
1100 FMUD=FLOAT (MUD) 

1110 THTU=THTDIF*CFMUD- 1 .0) 

1120 THTUP1 =THTDIF*FMUD 
1130 THTUPO=THTUP 1 -THTO 
1140 THTUO=THTU-THTO 

1 1 50 FINTTMCMUD, 1 , 1 ) = 1 .O/SGRTC FMUBAR) 

1160 D0 10 NTUD=2, NTUVMX 
1170 10 FINTTMCMUD, 1,NTUD)=0.0 
1180 D0 20 M TO D= 2, MTOVMX 
1190 FM TO D= F' L 0 A T C M TO D ) 

1200 AT0V=FMT0D-1 .0 

1210 ARGUP 1 =ATO V*THTUPO 

1220 ARGU=AT0V*THTU0 

1230 C=SQRTCFMUBAR/2.0)/CAT0V*PI ) 

1240 FVALUE=C*C SI NC ARGUP 1 )-SIN(ARGU) ) 

1250 IFC CABSCFVALUE) ) *LT. C 1 .OE-6) ) FVALUE=0.0 
1260 20 FINTTMCMUD, MTOD, 1 )=FVALUE 
1270 D0 40 MT0D=2, MTOVMX 
1280 FMTOD=FL0ATCMTOD) 

1290 ATOV = FMTOD- 1 .0 
1300 ARGUP 1=ATQV*THT UP 0 
1310 ARGU=ATOV*THTUO 
1320 D0 40 NTUD=2, NTUVMX 
1330 FNTUD=FL0ATCNTUD) 

1340 ATUV=C (FNTUD-1 • 0 ) *FMUEAR) /2 • 0 
1350 IF CATOV.EQ.ATUV) G0 T0 30 

1360 C=CAT0V/CAT0V**2-ATUV**2) ) /SQRTC CP 1**2) /FMUBAR) 
1370 AC= CFNTUD-1 «0)*PI 

1380 F VALUE=C* C COS C AC ) *SIN C ARGUP 1 ) - SI NC ARGU) ) 

1390 IFC CABSCFVALUE) ) «LT. C 1 .0E-6) ) FVALUE=0.0 
1400 FINTTMCMUD, MTOD, NTUD)=F VALUE 
1410 G0 T0 40 

1420 30 C=1 .O/SQRTC FMUBAR) 

1430 FVALU£ = C*C0S C ARGU) 

1440 IFC CABSCFVALUE) ) .LT. C 1 .OE-6) ) FVALUE=0.0 

1450 F I NTTM C MUD, MTOD,NTUD)=F VALUE 

1460 40 CONTINUE 

1470 1000 CONTINUE 

1480 WRITEC 1,2) F I NTTM, THTO 

1490 CALL CL0SEFC1) 
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TI T3 5 5 CONTINUED 


1500 PRINT*" THT0= "* THTO* "* i v il)BAR= "> MUBAR 

1510 p ft I NT* "MTQ Vi v iX = "* MTOVMX * "* NTUVMX = "* NTUVMX 

1520 STOP 
1530 END 
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READY 

0LD:TIT3bb 

HEADY 
000 1 $NDM 

RUN 


TIT35S 


19:46 


NAR JAN. 1 2> 1971 


INPUT: THT0*MUBAR,MT0VMX> NTUVMX 
?0 .00* 3* 5* b 

THT0= 0 .000000000E+00* MUEAR= 
MT0VMX= b* NTUVMX= b 

STOP 

RUNNING TIME: 2.3 SECS 1/0 TIME : 


D-33 


generat e T0 T}$$\ 


3 


.7 SECS 




1 


Tii 13 bb 


0000 1 
000 0 2 
00 00 3 
00004 
0000 3 
00006 
00 0 0 7 
00008 
0000 9 
000 1 0 
000 1 1 
000 1 2 
000 i 3 
000 1 4 
000 1 b 
000 1 6 
000 1 7 
000 18 
000 1 9 
00020 
00021 
00022 
0002 3 
000 2 4 
000 2 b 
000 26 


b • 7 73 50 2 7 E - 0 1 
3.3761862E-01 
1 . 6880931 E-0 1 
0 • 0 0 0 0 0 0 0 E + 0 0 
8 . 440 46 bbE-02 

0 • OOOOOOOE + QO 
3.81 97186E-01 
b . 4367409 E-0 1 
0 .OOOOOOOE+OO 
1 . 38898S6E-0 1 
0 .OOOOOOOE+OO 
b.9683104E-02 

1 « 9098 b9 3E-0 1 
b . 7 7 3 bO 2 7 E - 0 1 
2. 728370 b£-0 1 
0 .OOOOOOOE+OO 
2 • 4803368E-02 
5 « 8764902E-02 
0 .OOOOOOOE+OO 
4. 4937866E-0 1 
0 .OOOOOOOE+OO 
1 ♦ 3641 8b2E-02 
2.9841 552E-02 
0 .OOOOOOOE+OO 
9 . 5492966E-02 
0 .OOOOOOOE+OO 


b • 7 73 50 2 7 E~0 1 
-6 . 7 b 2 3 7 2 4 E - 0 1 
3* 3 7 61862 E- 01 
0 .OOOOOOOE+OO 
- 1 .6880931 £-01 
0 .OOOOOOOE + OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
1 . 1 936620E-01 
-3. 819 7 186 E- 03 
5. 773S027E-01 
- b • 456 7409 E-0 1 
0 -OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
0 .OOOOOOOE+OO 
2 • 7283705E-02 
-5.9683104E-02 
0 .OOOOOOOE+OO 
1 .9098593E-01 


5. 7 73 502 7E-0 1 
3.3761 8 6 2 £ - 0 1 
-1 -6880931 E-0 1 
0 .OOOOOOOE+OO 
8 . 4404655E-02 
0 -OOOOOOOE+OO 
-3. 8197186 E-01 

- 5. 4567409E-0 1 
0 .OOOOOOOE+OO 
1 • 3889886E-0 1 
0 .OOOOOOOE+OO 

-5.9683104E-02 
1 • 9098593E-0 1 
b. 7735027E-0 1 
2 • 7283 70 5E-0 1 
0 .OOOOOOOE+OO 
-2 . 4803368E-02 
5.8764902E-02 
0 .OOOOOOOE+OO 
-4.4937866E-01 
0 .OOOOOOOE+OO 

- 1 .3641 852E-02 
2.9841 552E-02 
0 .OOOOOOOE+OO 

-9 . 5492966 E -02 
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B3D353 


010GC THREE-DIMENSIONAL BAFFLED CHAMBER MAIN PROGRAM 
1000 DIMENSION F I NTRM ( 5* 3* 5* 3 > * F I NTTM ( 3* 5* 5) 

10 1 0& * ATOM ( 5 ) * AROM C 5* 3 ) * BESSOM (5*3).* ATUM ( 5) * ARUM ( 5* 3 ) * BESSUM C 5 * 3 ) 
1020& * HLOM ( 5*3)* HRUM (5*3) 

10 30& * ASOM( 5* 3) * VA0M( 5* 3)* BSUM(3* 5* 3) * VBUM ( 3* 5* 3) 

10 40& * BLOM (5*3) 

10 50 COMMON MUBAR* R 1 * ATOM* AROM* BESSOM* ATUM* ARUM* BESSUM 
1060& * N Rp NT* TVAL* NTPNT* RVAL 

10 70 COMPLEX BETAOP* BETA2L* El R1 ST* El R 1 SP* El R1 FD* El GNR 1 * El R 1 2 
1080& * AZOV* ARGO* RT ARGO* DELLOV* CTANO* AZOV* ARGU* RTARGU* DELRUV* CTANU 
10905: * HLOM* HLOV* HRUM * HRUV 

1 100& * ASOM* ASOV* VAOM* BSUM* BSUV* BLOM* BLOV* VBUM 
1 1 1 05 * TM* SERSUM* SUM 1 * SUM2* CHARFN* CHARFO* CHARFL* CHARFU 
1120 CALL 0PENF ( 1 * "RN5353”) 

1130 CALL 0PENF (2* "RT5353 1 ') 

1140 CALL 0PENF (3*"T0T355") 

11 50 2 FORMAT (4I5/I5*E10»2*I5) 

1160 6 FORMAT (F 5. 2* 3(/5X* 2F9 . 5) ) 

11 70 8 FORMAT ( IX* 5E1 2 • 5) 

1180 10 F0RMAT( IX* 3E1 5. 7) 

1190 11 F 0RMAT ( //"THREE-DIMENSIONAL BAFFLED CHAMBER”) 

1200 12 FORMAT ( //T2* "MUBAR= M * 12/ 

121Q&T6*"Ri=”,F5.2*T21*"Z0 = "*F5.2*T36*”Zl="*F5.2*T51*"Z2= e L F5.2/ 
1220&T1 1 *”BETA0P="* 2F 5 . 2* T3 1 * ”BETA2L= ”* 2F5.2) 

1230 14 FORMAT ( /T2* "TANGENTIAL MODE NUMBER"* 13) 

1240 16 FORMAT ( /T2* "MTO VMC = "* 1 2* T 1 6* "MRO VMC = "* I 2/ 

1250&T2* "NTUVMC = "* 12* T1 6* "NRUVMC="* I 2//T2* " JI T="* I 2) 

1260 18 FORMAT (//Tl 6* "El GNR 1"* T 52* "SUM 1 "/T 1 6 , "CHARFN"* T 52* "SUM2") 
1270 19 F0RMAT(T1O* 1PE16. 10*"*"* 1PE1 7. 10) 

1280 20 F 0RMAT (2( 1PE16. 10*"*"* 1 PE 1 6 • 10*" :"* 1 P E 1 6 • 10*"*"* 1 P E 1 6 • 10/) ) 
1290 READ ( 1*2) MTOVMX* MRO VMX* NTUVMX* NRUVMX* I DUM 1 * DUM 1 * I DUM2 
1300 READ ( 1*6) ( ATOM(MT) * ( AROM ( MT* MR) * BESSOM ( MT* MR) * MR= 1 * MROVMX ) * 
1310&MT=1 *MTOVMX) 

1320 READ ( 1*6) (ATUM (NT)* (ARUM ( NT* NR) * BESSUM ( NT* NR) * NR= 1 * NRUVMX > * 
1330&NT=1 * NTUVMX) 

1340 READ( 2*8) FINTRM 

1350 READ( 3*10) F I NTTM* PHEO OP 

1360 1001 PRINT*" 1 INPUT : MUBAR* R1*Z0*Z1*Z2*BETA0P* BETA2L" 

1370 INPUT* MUBAR* R1 *Z0*Z1* Z2* BETAOP* BETA2L 

1380 1002 CONTINUE 

1390 NM0DE= 1 

1400 BSUMO 1 = + 0 • 5 

1410 BSUMO 2 = - 1 .0 

1420 B S UM 0 3 = + 0 • 5 

1430 1003 PRINT*" 3 I NPUT :MT0VMC* MROVMC* NTUVMC* NRUVMC* JI T" 

1440 I NPUT* MTO VMC* MRO VMC* NTUVMC* NRUVMC* JI T 
1450 1004 PRINT*" 4 INPUTsITYPE" 

1460 INPUT* ITYPE 
1470 I DETAL=0 

1480 IF(ITYPE) 1008*1005*1006 


B-35 



2 


B o T 3 3> ! i 0 > \ T I 3; U Z D 


1490 100 5 PRINT;" 5 INPUT: EIR 1 ST; El R 1 SP; NEI R1 S> I DETAL" 
1 b 0 U 1 N P U T > E I R 1 ST ; El R 1 SP; NEI R 1 S; I DETAL 
1 b 1 0 IF ( I DETAL • NE • 1 ) GO T0 100 7 

1 5 2 0 P H I N 7 ; ’ ’ INP U T :IPBSUM;I P VBUM; I PASO M ; IPVAOM" 

1 0 30 I NP UT ; I P BS UPl > IP VBUM; IPASOM; IPVAOM 
I 040 f ;< I N 7 > *' I NPSJI MP8SUM; MP VBUM; MPASOM; MPVAOM" 

1 5b| I NPU T; MPBSUM;MPVeUM;MPASOM;MPVAQM 

1 b 6 0 G 0 T 0 100 7 

1570 1006 PRINT;" 6 I NP UT : E I R 1 ST; El R 1 FD; NEI R 1 M, CHFNTL" 
1580 INPUT; El SI ST; El R 1 FD; NEI R 1 M; CHFNTL 

1590 1007 CONTINUE 
1600 PRINT 11 

1610 PRINT 1 2; MUBAR; R1 ; ZO; Z 1 ; Z2; BETAOP; BETA2L 

1 6 2 0 p R I N I' 1 4 j NM 0 D E 

1630 -HINT 1 6; MTOVMC; MROVMC; NTUVMC; NRUVMC; JI T 

i 6 40 P R I N 3 1 8 


17G0C .DEFI.NTITI0NS 
1 ? I 0 C 
1720 ,i| « ; 0: 

1730 EIR12=21GNR1**2 
1740 D0 30 MT0D=1;MT0VMC 
1 7 00 D 0 2 0 M RO D= 1 ; M RO VM C 

1760 IF ( (MTOD.EQ. 1 ) .AND. CMROD.EQ. 1 ) ) G0 T0 29 

1 7 7 0 A K 0 V = A R 0 M ( M T 0 D ; M RO D ) 

1 7 8 0 A R G V 2 =ARQV * * 2 
1 7 9 0 R T A R GO = E I R 1 2 - A RO V 2 
1 8 0 0 A Z 0 V = C3QRT ( RTARGO) 

18 1 0 AkGO=AZOV*CZ 1 -ZO ) /R1 

1820 DELLOV=CO .0; 1 .0>*EI GNR 1 *BETAOP/AZO V 

18 30 C T A N 0 = C S I N ( A R G 0 ) / C C 0 S C A R G 0 ) 

1840 HLOV=-AZOV*( CTANQ+DELLOV ) / ( 1 .0-DELL0V*CTAN0> 

1 8 5 0 B L 0 Pi C M T 0 D ; M R 0 D ) = H L 0 V / R 1 

1860 29 CONTINUE 

1870 30 CONTINUE 

1880 HLQMC 1 ; 1 )=0 .0 
1890 D0 40 NTUD= 1 ; NTUVMC 

1900 DS 39 NRUD= 1 ; NRUVMC 

19 10 AfiUV=ARUM C NTUD; NRUD) 

1 9 2 0 A R U V 2 = A R U V * * 2 

1930 RTARGU= El R 1 2- ARUV2 
19 40 AZ U V = CSO RT CRTARGU) 

1950 ARGU=AZUV*CZ2-Z1 )/Rl 

1960 DELRUV= (0*0; 1 . 0 ) *EI GNR 1 *BETA2L/AZUV 

1970 CTANU=CSIN(ARGU) /CC0SCARGU) 

1980 HRUV=+AZUV*<CTANU+DELRUV)/C 1 .0-DELRUV*CTANU> 
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19 9 0 H R U M CNTUL)*N H U 0 ) = H R UV / R 1 
2000 39 CONTINUE 

2010 40 CONTINUE 
2020 DO 43 MUD= 1 >MUBAR 
2030 DO 42 NTUD= 1 > NTU9MC 
2040 DO 41 N.RUD= 1 > NRUVMC 
2050 6SUMCMUD* NTUD* NRUD)=0 -0 
2060 41 CONTINUE 

20 70 42 CONTINUE 
2080 43 CONTINUE 

2090 BSUMC 1 , 1 , 1 )=BSUM01 
2100 BSUM ( 2 * 1 * 1 )= BSUM 02 
2110 BSUM ( 3.» 1 .» 1 ) =BSUMQ3 
2120 45 CONTINUE 

21 30C 

2140C BLOM CALCULATION 

21 50C 

2160 DO 90 M TO D= 1 > M TO VM C 
2170 DO 80 M R 0 D = 1 * M RO V M C 

2180 IFC CMTOD.EQ. 1 ) .AND. CMROD.EQ. 1 ) ) GO TO 80 

2190 SERSUM=0 • 0 

2200 D0 70 MUD = 1 > MUBAR 

2210 DO 60 NTUD= 1 > NTUVMC 

2220 DO 50 NHUD= 1 > NRUVMC 

2230 TM=BSUMC MUD> NTUDj NRUD) *HRUM C NTUDj NRUD) 

2240 VBUM C MUD j NTUD j NRUD)=TM 

22 50 TM=TM*FI NTRMC MTODj MROD.* NTUDj NRUD) *F I NTTM (MUD.* MTO D* NTUD) 
2260 S E R S UM = S E R S UM + TM 

2270 50 CONTINUE 

2280 60 CONTINUE 

2290 70 CONTINUE 

2300 BLOM C MTO Dj MROD) =SERSUM 

23 1 0 VAO M C MTO Dj MROD) = SERSUM 

2320 ASOM C MTO D.* MROD ) =SERSUM/HLOM ( MTOD* MROD) 

2330 80 CONTINUE 
2340 90 CONTINUE 

23 50 BLOMC 1 * 1 )=0 -0 
2360 VAOMC 1 , 1 )=0 *0 
23 70 ASOM C 1 ., 1 )=0 .0 

2380 IFCJI .GE.JIT) GO TO 150 

2390 JI=JI+1 

2400C 

241 OC ITERATION EQUATION CALCULATION 
2420C 

2430 DO 140 MUD= 1 > MUBAR 
2440 D0 130 NTUD= 1 > NTUVMC 
2450 DO 120 NRUD= 1 > NRUVMC 
2460 SERSUM = 0 « 0 
2470 DO 110 MTOD= 1 > MTOVMC 
2480 DO 100 MROD= 1 > MRQVMC 
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3 CONTI NUED 


2490 IF C ( MTO D « EG « 1 ) . AND • ( MROD • EG) « 1 ) ) GO T0 100 

2509 TM = As Jrt( MTOD* MROD) 

2 1 0 f M = T M *FINTR M C M T 0 D * M R 0 D* NTUD* NRUD) *F I NTT M ( M U D , M T 0 D , N I U D ) 

2 . 2 0 S ER S UM = 5 E A S UM + 1 A 

c',b3u i DU CON TINUE 

2d 40 1 i 0 C b N T I N U E 

2ddO &3UM ( MUD* NTUD* NRUD) = SERSUM 

2560 120 CONTINUE 

2 5 7 0 130 C 0 N TINUE 

2580 140 CONTINUE 

2590 GO TO 45 

2600 150 CONTINUE 

26 1 OC 

2620C CHARACTERISTIC EQUATION CALCULATION 

26 3 0 C 

26 40 S E R S U M = 0 • 0 

2650 DO 170 MTO D= 1 * MTOVMC 

2660 DO 1 60 MR0D= 1 > MROVMC 

2670 IRC ( NTOD.EQ • 1 ) .AND. CMROD.EQ. 1) ) G0 T0 160 

2680 ASO 7=AS0MCMT0D*MR0D) 

2690 TM=C AS0V**2)*HL0M(MT0D*MR0D) 

2 • 0 0 S E R S UM = S E R S U M +TM 

2710 160 CONTINUE 

2720 170 CONTINUE 

2/30 SUM 1 =SERSUM 

2740 SEES UM=0 » 0 

2750 DO 200 MUD= 1 * MUBAR 

27 60 DO 1 90 NTUD= 1 > NTUVMC 

2770 DO 180 NRUD= 1 * NRUVMC 

2780 BSUV=BSUM(MUD* NTUD* NRUD) 

279 0 TM=< BSUV**2)*HRUMCNTUD* NRUD) 

28 0 0 S E R S UM = S E R S UM + TM 
2810 ISO CONTINUE 
2220 190 CONTINUE 
2830 200 CONTINUE 
28 4 0 S U M 2 = SERSU M 
28 50 CHARF N = SUM 1 - SUM 2 

28 6 0 C 

2b 7 0 C G £ N E R A L PRINTOUT 
2880 C 

£890 IF ( NCNTRL » EQ • -2) PRINT 20* El GNR 1 * SUM 1 * CHARFN* SUM 2 

29 0 0 I F ( I DETAL . NE . 1 ) GO TO 260 

29 1 0 C 

2920C DETAILED PRI NT0UT- - BSUM* VBUM* ASOM* VAOM 

29 30 G 

2940 IF < C IPBSUM+MPBSUM) .NE.O) PRINT* t2*"BSUM"*t 

2950 IFC IPBSUM.NE. 1 ) GO TO 241 

29 60 PRINT I 9* C C CBSUMCMD* NT* NR) *MD= 1 * MUBAR) * NT= 1 * NTUVMC )* NR= 1 * NRUVMC) 

2970 241 CONTINUE 

2980 IF CMP 3 SUM . EQ . 1 ) CALL PV3DBCC BSUM* NTUVMC* NRUVMC) 


D-38 



5 


B3D353 CONTINUED 


2990 IFC C IPVBUM+MPVBUM) .NE.O) PRINT# »2# "VBUM"# t 
3000 IFCIPVBUM.NE. 1 ) G0 T0 243 

3010 PRINT 1 9# C ( C VBUMCMD# NT# NR) #MD= 1 # MUBAR) # NT= 1 # NTUVMC) # NR= 1 # NRUVMC) 
3020 243 CONTINUE 

3030 IFCMPVBUM.EQ.l ) CALL PV3DBCC VBUM# NTUVMC# NRUVMC) 

3040 IF C C IPAS0M+MPA50M) .NE.O) PRINT# t2#"AS0M"# t 
3050 IF C IPASOM . NE • 1 ) G0 T0 245 

3060 PRINT 19# C ( ASOM ( MT# MR ) # MT= 1 # MTO VMC ) # MR= 1 # M RO VM C ) 

3070 245 CONTINUE 

3080 IFCMPASOM.EQ. 1 ) CALL PV3DMCC ASOM# MTO VMC# MROVMC# PHEO 0P ) 

3090 IFCC IPVAOM+MPVAOM) .NE.O) PRI NT# t 2# " VAOM"# t 
3100 IFC IP VAOM • NE • 1 ) G0 T0 247 

31 10 PRINT 19# C CVAOMCMT#MR) # MT= l # MTO VMC) #MR=1 #MROVMC) 

3120 247 CONTINUE 

3130 IFC MP VAOM . EQ • 1 ) CALL PV3DMCC VAOM# MTOVMC# MROVMC# PHEO0P ) 

3140 260 CONTINUE 
31 50 C 

3160C CHANGE El GNR 1 - - STEP CHANGE 0R R00T FINDER 
31 70C 

3180 NEI R 1 D=NEI R1 D+ 1 

3190 IFCITYPE-EQ. 1) G0 T0 300 

3200 EIGNR1=EIGNR1+EIR1SP 

3210 IFCNEIR1D.LT. NEIR1S) G0 T0 1 

3220 G0 T0 1008 

3230 300 CONTINUE 

3240 El R 1 FR=E I R 1 FD 

3250 EIR1FI=-C0.0# 1 .0)*EIR1FD 

3260 IF CABS CEIR1FR) .GT.C1 .0E-10)) G0 T0 301 

3270 CHARFI =-C0*0# 1 *0) *CHARFN 

3280 CHARFN=+CO .0# 1 .0)*CHARFI 

3290 G0 T0 302 

3300 301 CONTINUE 

3310 IFCABSCEIR1FI) .GT.C1 .0E-10)) G0 T0 302 

3320 CHARFR=CHARFN 

3330 CHARFN=CHARFR 

3340 302 CONTINUE 

3350 NCNTRL=NCNTRL+ 1 

3360 IFC NCNTRL ) 310#320#330 

3370 310 CHARFO=CHARFN 

3380 IFCCCABSCCHARFO) .LT.CHFNTL) .0R.CNEIR1D.GT.NEIR1M)) G0 TO 1008 

3390 El GNR1 = EI GNR 1 - El R1 F D/2 • 0 

3400 GO T0 1 

3410 320 ,CHARFL=CHARFN 

3420 EIGNR1=EI GNR 1+EIR1FD 

3430 GO TO 1 

3440 330 CHARFU=CHARFN 

3450 El GNR 1 = El GNR 1 - El R1 FD/2 «0-CHARF0*EI R1 FD/C CHARFU- CHARFL) 

3460 NCNTRL=- 2 

3470 GO TO 1 

3480 1008 PRINT# »2# »2 
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B3 D 5 53 CONTINUED 


349 OC 

3500 C WHERE TO G£? 

33 10 0 

3540 PRINT, " INPUT: I GET©” 

3530 INPUT, I GOT© 

3540 GO TO C 1001, 1002, 1003, 1004, 1005, 1006,9999), I GOTO 

3530 9999 STOP 
3560 END 
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PV3DBC 


5000 SUBROUTINE PV3DBCC C0EF* NTUVMC* NRUVMC) 

5010 DIMENSION C0EFC3* 5* 3) > REFERM (5,3), BESRVM C 5* 3) *PVVALC 50) * RTVAL (. 50 ) 
5020& > ATOM ( 5) * ARQM ( 5* 3) > BESSOM C 5* 3) * ATUM ( 5) > ARUMC 5* 3)> BESSUMC 5* 3) 

50 30& * COSFNC 5) *MUDM( 50 ) 

50 40 COMMON MUBAR* R1 > ATOM# AROM* BESSOM* ATUM* ARUM* BE S S UM 
50 50 & * NRPNT* TVAL* NTPNT* RVAL 
50 60 COMPLEX C0EF 

5070 5 FORMAT ( /5X* "R- VARIATION COMPARTMENT'S I 3* 5X* "TVAL=’S F 1 0 • 4* iC *P I "/) 
5080 10 F0RMATC2X* 1P2E1 7. 10) 

5090 15 F0RMATC/5X* "THETA-VARIATION'S 5X*"RVAL = 'S 1PE1 7. 10/) 

5100 20 F0RMATC2X* 1PE1 7. 1Q*QPF10 • 4* "*PI’S I 10) 

51 10 PI=3 • 1 41 5926536 

5120 PIT2=2.0*PI 

5130 FMUBAR = FLOAT CMUBAR) 

5140 P I 0FMU=P I /FMUBAR 
5150 PI T20F = P I T2/FMUBAR 
5160 R1202=CR1**2)/2.O 
5170 DO 110 NTUD= 1 * NTUVMC 
5180 DO 100 NRUD= 1 > NRUVMC 
5190 EP SV=P 1 0FMU 

5200 IFCNTUD.EQ. 1 ) EPSV=PIT20F 
5210 ATUV2= ( ATUM ( NTUD) ) 

5220 ARUV2=CARUMCNTUD*NRUD) )**2 

5230 REFRV2=EPSV*R1202*( 1 • 0- C ATUV2/ARUV2) ) 

5240 8ESUV2=CBESSUMCNTUD*NRUD) )**2 

5250 REFERMCNTUD*NRUD)=SQRTCREFRV2*BESUV2) 

5260 100 CONTINUE 
5270 110 CONTINUE 

5280 REFERM Cl* 1 ) = SORT (P I T20F*R1 202) 

5290 1 PRINT* T2*t," INPUT: NRPNT * TVAL* NTPNT* RVAL * T STOP'S t * 

5300 INPUT* NRPNT* TVAL* NTPNT * RVAL* TST0P 

5310 IFC(NRPNT.EQ.O) .AND* (NTPNT. EQ.O) ) RETURN 

5320 IF ( NRPNT .EQ.O) G0 T0 1000 

5330 TVAL0P=TVAL 

5340 TVAL=TVAL*PI 

5350 IFC (TVAL. GT. 0.0) .AND. CTVAL.LT. PIT20F) ) MUD=1 

5360 IFC (TVAL . GT.PIT20F) . AND . ( TVAL .LT • ( 2 . 0*P I T20F ) ) ) MUD-2 

53 70 IFC CTVAL.GT. (2.O+PIT20F) ) .AND* (TVAL .LT»(3»O*PIT20F) ) ) MUD=3 

5380 FMUD=FL0ATCMUD) 

5390 D0 200 NTUD= 1 * NTUVMC 

5400 C0SFNCNTUD)=C0SCATUMCNTUD)*CTVAL-CFMUD-1 .0)*PIT2BF) ) 

5410 200 CONTINUE 

5420 RI'NC = R1 /NRPNT 

5430 NRPNT 1 =NRPNT+ 1 

5440 RP0S=O.O 

5450 D0 500 NRD= 1 * NRPNT1 

5460 RTVAL ( NRD) =RP0S 

5470 SUM = 0 «0 

5480 D0 410 NTUD=1* NTUVMC 
5490 D0 400 NRUD=1* NRUVMC 
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PVo DSC CON 1 I NUED 


5500 RPUSER=ARUM( NTUD* NRUD) *RP0S/R1 

55 10 CALL B E S F Np CAT UM (NTUD) * RPBSER* BESSFN) 

5520 SU; v i = SUM + C0£F (MUD* NTUD> NRUD) *C0SFNC NTUD) *BESSFN/REFERMC NTUD* NRUD) 

5530 400 CONTINUE 
5 40 410 CONTINUE 
5 5 50 P 0 V A L C N R D ) = S UM 
5560 RP OS - RP 03 + R I NC 
5570 dOO CONTINUE 
5560 PRINT 5* MUD* TVAL0P 

5590 PRINT 1 0 * (PWAL ( NRD) * RTVAL ( NRD) * NRD= 1 * NRPNT 1 ) 

5600 1000 CONTINUE 
5610 IF ( NTPNT . EQ • 0 ) G0 T0 1 
5620 DO 1110 NTUD = 1 * NTUVMC 
5630 DO 1100 NRUD=1 *NRUVMC 

56 40 RP 0SER=ARUM C NTUD* NRUD) *RVAL/R1 

5650 CALL BESFNP C ATUM C NTUD) * RP0SER* BESSFN) 

5660 BE S RV'M ( NTUD* NRUD) =BESSFN/REFERMC NTUD* NRUD) 

56 70 1 100 CONTINUE 
5660 1110 CONTINUE 
5690 II NC=P I T2/NTPNT 
5700 TI NC0P=2 .O/NTPNT 
5710 NTPNT I =NTP NT +1 
5720 TP0S= 1 .OE- 10 
5730 TP0S0P=O.O 
5740 TST0P=TP0S+TST0P*PI 
5750 D0 1 500 NTD= 1 * NTPNT1 
5760 RTVALCNTD) =TP0S0P 

5770 IFCCTP0S. GT. 0.0) .AND. CTP0S.LT. PIT20F)) MUD=1 

5780 IFC CTP0S .GT.PIT20F) * AND » ( TP0S «LT « C 2 • 0*PI T20F ) ) ) MUD=2 

5790 IFC C TP0S • GT . C 2 • 0*P I T20F ) ) • AND • C TP0S *LT • C 3 . 0*P I T20F ) ) ) MUD=3 

5800 M U DM C N T D ) = M U D 
5810 FMUD=FL0ATCMUD) 

5820 SUM=0.0 

58 30 DO 1410 NTUD= 1 * NTUVMC 
58 40 DO 1 400 NRUD= 1 * NRUVMC 

5850 TM=BESRVM( NTUD* NRUD) * COS (ATUM C NTUD) * C TP0S~ C FMUD- 1 . 0 ) *P I T20F ) ) 
5860 SUM=SUM+TM*C0EF( MUD* NTUD* NRUD) 

58 70 1400 CONTINUE 
5880 1410 CONTINUE 

58 9 0 PVVAL C NTD) =SUM 
5900 NTMAX=NTD 

5910 IFCTP0S.GT.TST0P) G0 T0 1600 

5920 TP05=TP0S+TINC 
5930 TP0S0P=TP0S0P+TI NC0P 

59 40 1500 CONTINUE 
5950 1600 CONTINUE 
5960 PRINT 1 5* RVAL 

5970 PRINT 20* ( P VVALC NTD) * RTVAL C NTD) * MU DM ( NTD > * NT D= 1 * NTMAX ) 

5980 G0 T0 1 
5990 END 
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6000 SUBROUTINE P V3DMC C C0EF, M TO VMC, MROVMC, PHEO 0P ) 

60 1 0 DI MENS I ON C0EFC5,3), REFERM C5,3), BESRVM (5,3). PVVAL C 50 ) , RTVALC 50 ) 
6020& , ATOM C 5) , AROM (5,3), BESSOMC 5, 3) , ATUMC 5) , ARUM ( 5,3), BBS SUM C 5, 3) 

60 3 0 & , C0SFNC5) 

6040 COMMON MUBAR, R 1 , ATOM, A ROM, BESSOM, ATUM, ARUM, BESSUM 
6050& , NRPNT, TVAL, NTPNT , RVAL 
6060 COMPLEX C0EF 

60 70 5 F0RMATC/5X,"R- VARIATION MAIN CHAMBER", 7X, "TVAL= ", FI 0 . 4, "*P I’ '/) 
6080 10 F0RMATC2X, 1P2E1 7. 10) 

6090 1 5 F0RMATC/5X, "THETA-VARIATION", 5X,"RVAL = ", 1PE1 7. 10/) 

6100 20 F0RMATC2X, 1PE1 7. 10,0PF10 . 4, "*PI") 

61 10 P I = 3 • 1 41 5926536 
6120 PIT2=2.0*PI 

6130 R1202=CR1**2)/2.O 
6140 DO 110 MTO D= 1 , MTOVMC 
6150 DO 100 MROD= 1 , MROVMC 
6160 EPSV=P I 

6170 IFCMTOD.EQ. 1 ) EPSV=PIT2 

61 SO AT0V2=CAT0MCMT0D))**2 

6190 AR0V2= < AR0MCMT0D,MR0D) )**2 

6200 REFRV2=EPSV*R1 202* ( 1 .0- C AT0V2/AR0V2) ) 

6210 BESO V2= ( BESSOM (MTOD, MROD) ) **2 

6220 REFERMC MTO D,MR0D)= SORT CREFRV2*BES0V2) 

6230 100 CONTINUE 
6240 110 CONTINUE 

6250 REFERMC 1, 1 ) =SQRTCPI T2*R1 202) 

6260 PHEEO=PHEO0P*PI 

6270 1 PRINT, »2,t," I NPUT : NRPNT, TVAL, NTPNT, RVAL, TST0P*', t * 

6280 I NPUT, NRPNT, TVAL, NTPNT, RVAL, TST0P 

6290 IF ( ( NRPNT • EQ ♦ 0 ) .AND • ( NTPNT • EQ • 0 ) ) RETURN 

6300 IF (NRPNT -EQ «0) G0 T0 1000 

6310 TVAL0P=TVAL 

6320 TVAL=TVAL*P I 

6330 D0 200 MT0D= 1 , MTOVMC 

6340 C0SFNC MTOD) = COS ( ATOM CM TOD) * ( TVAL-PHEEO ) ) 

6350 200 CONTINUE 

6360 RI NC=R 1 /NRPNT 

6370 NRPNT 1 =NRPNT+ 1 

6380 RP0S=O *0 

6390 DO 500 NRD= 1 , NRPNT 1 

6400 RTVALC NRD)=RP0S 

6410 SUM = 0 • 0 

6420 DO 410 MT0D=1, MTOVMC 

6430 D0 400 MROD= 1 , MROVMC 

6440 RP0SER=AROMCMTOD,MROD)*RP0S/R1 

6450 CALL BESFNP C ATOM (MTOD) , RP0SER, BESSFN) 

6460 SUM=SUM + C0EF ( MTOD, MROD) *C0SFN(MTOD) *BESSFN/REFERM C MTOD, MRO D) 

6470 400 CONTINUE 
6480 410 CONTINUE 
6490 PVVALC NRD) =SUM 
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PV 3 iJhiC CONTINUED 


6 5 u 0 H p fc 5 - R P b S + R I N C 
6 1 50 0 CONTINUE 

6 5 2 5 P H I ! v I TVAL0P 

6530 PRINT 1 Oj (PV/VALCNRD) .» RIVAL ( NRD) , NRD= 1 * NRPNT1 > 

6b 40 1000 CONTINUE 

6 5 50 I F C NTP N T * EQ • 0 ) G0 T0 1 

6560 D0 1110 MTOD= 1 * MTOVMC 

6570 DO 1100 MROD= 1 jMROVMC 

6580 R P 0 3 E R = A R 0 M ( M T 0 D > M R 0 D > * R V A L / R 1 

6 590 CALL BESFNPC ATOM (MTOD) * RP0SER> BESSFN) 

6 6 0 0 B E S H V M ( MTO Dj MROD) = BESSFN/REFERM C MTODj MRO D) 

6610 1 1 00 CONTINUE 
6620 1110 CONTINUE 
6630 TI NC=PI T2/NTPNT 
6640 TI NC0P=2 .O/NTPNT 
6650 NTPNT1 =NTPNT+1 
66 60 TP 0 S = 1 » 0 E ~ 1 0 

66 70 TP 3 SOP = 0.0 

6680 TSTOP=TP0S+TST0P*PI 
669 0 DO 1500 NTD= 1 * NTPNT1 
6700 RTVAL( NID) =TP0S0P 
6710 SUM = 0 « 0 

6720 DO 1410 MTOD= 1 j MTOVMC 

6730 DO 1400 MROD= 1 > MROVMC 

6740 TM= 3ESRVM C MTO D* MROD) *C0SC ATOM C MTOD) *CTP0S-PHEEO) > 

67 50 S UM = SUM + TM*C0EF CMTODjMROD) 

6760 1400 CONTINUE 
6770 1410 CONTINUE 
6780 PVVALCNTD) =SUM 
6790 NTMAX=NTD 

6800 IFCTP0S.GT.TST0P) G0 T0 1600 

6810 TP05=TP0S+TINC 
6820 TP 0S0P = TP0S0P + TI NC0P 

6830 1500 CONTINUE 
6840 1600 CONTINUE 

68 50 PRINT 1 5* RVAL 

6860 PRINT 20 j ( PVVAL ( NTD) * RTVAL (NTD) * NTD= 1 > NTMAX ) 

6870 G0 T0 1 
6880 END 
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7000 SUBROUTINE BESFNP ( ORDER* ARGMNT* BESFNV) 
7010 I 0RDER= IFIX(ORDER) 

7020 0RDDIF=0RDER-FL0AT(I ORDER) 

7030 IFC0RDDIF . GT. C 1 .0E-8) ) G0 T0 1000 
70 40 C INTEGER ORDER BESSEL FUNCTIONS 
7050 IFCI0RDER.EQ-O) BESFNV=1.0 
7060 IFC I0RDER.NE-O) BESFNV=0.0 
7070 IFCARGMNT.LT. ( 1 .OE-8) ) RETURN 
7080 CALL BESSOPCARGMNT*BESSOV) 

7090 IFCI0RDER.GT.O) G0 T0 100 

7100 BESFNV=BESS0V 

7110 RETURN 

7120 100 CONTINUE 

7130 CALL BESS 1 P C ARGMNT * BESS 1 V) 

7140 IF C I ORDER • GT • 1 > GO T0 200 

7150 BESFNV=BESS 1 V 

7160 RETURN 

7170 200 CONTINUE 

7180 BESL 1 =B£SS0V 

7190 BES0 = BESS 1 V 

7200 I 0RD= 1 

7210 300 CONTINUE 

7220 0RDT2=2.O*FL0ATCI0RD> 

7230 BESF NV= C 0RDT2/ARGMNT) *BESO-BESL 1 

7240 I 0RD= I 0RD+ 1 

7250 IFCI0RD.EQ.I ORDER) RETURN 

7260 BESL1=BES0 

7270 BESO=BESFNV 

7280 60. T0 300 

7290C HALF- INTEGER BESSEL FUNCTIONS 
7300 1000 CONTINUE 
7310 PI=3. 1 41 5926536 
7320 BESFNV=0»0 

7330 IFCARGMNT.LT. C 1 .0E-8) ) RETURN 
7340 CHALF= SORT C2. 0/CPI *ARGMNT> ) 

7350 BESM1 2=C0SC ARGMNT) 

7360 BESP12=SINCARGMNT) 

7370 IFCI0RDER.GT.O) G0 T0 1100 

7380 BESFNV=CHALF*BESP12 

7390 RETURN 

7400 1100 CONTINUE 

7410 BESL 1 =BESM 1 2 

7420 BES0=BESP 1 2 

7430 I 0RD= 1 

7440 1200 CONTINUE 

7450 0RDT2=2.O*CFL0ATCI0RD)-O.5O> 

7460 BESFNV= C 0RDT2/ARGMNT) *BES0~BESL 1 
7470 IFCI0RD.NE* I ORDER) G0 T0 1300 
7480 BESFNV=CHALF*BESFNV 
7490 RETURN 
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bESFNP CONTINUED 


750 0 

1300 CONTINUE 

/ D 1 0 

IC'KD = I0RD+1 

75 20 

b E 3 L 1 = b E 3 U 

7530 

BE1S0 = BESF NV 

7 5 A 0 

GO T 0 1200 

7 3 30 

£00 
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BE' 50ft 


60 0 0 

3 U B 

KL-UT I NF BES so 

ft ( A K Ci 

60 1 0 

D I M 

EivS I 'itv CbO C 7 ) 

; CFO C 

60 20 

DA T 

-i O 0/1 • 0 y “ 2 • 

249 9 9’ 

60 30 4; 

y — 

3 • '9 4 4 4 E - 3 » 8 • 1 

E - 4 / , 

60 404 

j — 

9 • 8 1 2 E - 3 j 1 • 3 7 

23 7E- 

oO 50 & 

j - 

4 • 1 6 639 / E — 8 j — 

3 • 9 5 4 

80 6 0 4 

> 1 

• 3 5 3 8 E ~ 4 / 


80 7 0 

IF C 

A n G M N I • G f • 3 • 0 ) G G 

8 080 C 

A ft 

i J !«’ S T l < 6, 

N 0F< 

60 9 0 

A — A 

*■■■ 3 1 / > i\' i / 3 • 0 


8100 

SUN 

o 

• 

o 

it 


6 1 1 0 

u 

100 J=l»7 


8 1 2 0 

K= 2 

4 C J - 1 ) 


8130 

SUM 

= SUi-i + CBO ( J ) * ( 

X 4 4 K ) 

61 4u 

1 0 0 

C 0. N T I N U E 


c 1 50 

£ E 3 

S0v'= SUM 


6 1 6 0 

RETURN 


6 i / 0 C 

/O i-^ 

OMNI GREATER 

THAN 


T) , CTO ( 8 ) 

S' 7 j 1 • 2 6 06 20 b j 
CFO / 7 • 9 708456 
3j - 7 • 26 0 5 F. ~ 

E - 5 * 2 # 6 2 5 7’ 3 £ - 

TO 1000 
EQUAL TO 3*0 


- ■ j • ,j 

£- 1 ^ 
1 • 44 
3 } - 5 


1 6 38 

- 7 • 7 
7 6 E - 

• 4 1 ? 


t. o j 4 • 

E - 76 - 

4/6 CT 
5E- 4 j 


;■> • 
.V 


4 /V 

5 2 7 4 * 

. 5 • L j ' 

* 9 3 3 ' 


o • 


o 


6 1 8 0 1000 C 0 N T I N U E 
8190 a = 3 • 0 / A ft GW N T 
8200 SuV) = 0 *0 
8210 DO 1100 J=h7 
8220 K= J - 1 

6230 S U i v j — SUM + CFO ( J ) -! ; ( X 'TH-'K ) 

8240 1100 CONTINUE 
82 80 F0= SUM 
6260 SUrt = 0 »0 
82 7 U DO 1200 J=1>B 
826 0 K = J- 1 

8290 S U M =SUM + CT0< J>*< X * * K ) 

830 0 1200 CONTI NU E 
6310 T0 = SUM/X 

6320 bESSOV- C FO^COS CTO) ) /SORT ( ARGMNT) 

6330 ft F. TURN 
6340 END 
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BE S SIP 


9000 SUBROUTINE BESS IP < ARGMNT* 8ESS1 V> 

9010 DIMENSION CB 1 C 7) * CF 1 C 7) * CT 1 C 8 ) 

9020 DATA CB 1 /5 • 0 E- 1 * - 5 . 624998 5E- 1 * 2 • 1 0935 73E- 1 * - 3 . 9 54289 E- 2 
90 30& * 4. 43 31 9E-3*-3. 1 761 E-4* 1 . 1 09E- 5/* CF 1 / 7. 9 788 456E- 1 * 1 . 56 E- 6 
90 40 & * 1 . 659667E-2* 1 . 71 05£- 4* -2. 49 51 IE- 3* 1 . 1 3653E- 3* -2 • 0033E- 4 
90 SO & /*CT1 /3 .0* -2*3561 9449, 1 . 24996 12E- 1 * 5.6 5E-5* - 6. 3 78 79 E- 3 
9060& , 7.4348E-4* 7 . 9S24E- 4* -2 . 9 1 66E- 4/ 

9070 IF C ARGMNT. GT. 3.0) GO T0 1000 

9080 C ARGMNT LESS THAN 0R EQUAL T0 3.0 

9090 X = AR GMNT/3 » 0 

9100 SUM = 0 • 0 

9110 D0 100 0=1.. 7 

9120 K=2*(J~1> 

9130 SUM=SUM+CB1 ( J)*CX**K) 

9140 100 CONTINUE 
9150 EESS1 V=ARGMNT*SUM 
9160 RETURN 

91 70 C ARGMNT GREATER THAN 3-0 

9180 1000 CONTINUE 

9190 X=3 .0 /ARGMNT 

9200 SUM =0.0 

9210 DO 1100 J=1.7 

9220 K = J - 1 

9230 SUM=SUM+CF1 C J)*CX**K> 

9240 1100 CONTINUE 
9250 F 1 =SUM 
9260 SUM=0 • 0 
92 70 D0 1200 J=H8 
9280 K=J-1 

9290 SUM=SUM+CT1 ( J)*<X**K> 

9300 1200 CONTINUE 
9310 T1=SUM/X 

9320 BESS1V=(F1*C0SCT1 ) ) /SQRT (ARGMNT) 

9330 RETURN 
9340 END 
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1 INPUT:MUBAR,R1,Z0,Z1,Z2,B£TA0P,B£TA2L 
? 3, 1 . 0/0. 0*2.1 >2. 8,0. 0,0. 0,0. 0*0.0 

3 input:mtovmc,mrovmc,ntuvmc, nruvmc, jit 

3) 10 

4 INPUT: I TYPE 
?0 

5 INPUT:EIR1ST,EIR1SP,NEIR1S,I DETAL 
?1 . 10, 0*00*0 . 10, 0 .00, 5,0 


THREE-DIMENSIONAL BAFFLED CHAMBER 


MUBAR= 3 

Rl= 1.00 Z0 = 0.00 Z 1 - 2-10 Z2= 2. BO 

BETA0P= 0.00 0.00 BETA2L= 0.00 0.00 

TANGENTIAL MODE NUMBER 1 

MTO VMC= 5 MROVMC= 3 

NTUVMC= 5 NRUVMC= 3 

JIT=10 


EIGNR1 SUM1 

CHARFN SUM 2 

1 . 10000000QE+00, 0 .OOOOOOOOOE+OO : 2.9 43 4776 60 E- 03, 6 . 1 462900 1 6E- 1 5 
6. 13866 5321 E-0 3, 1 .40 76501 9 7E- 14 : -3 . 1 9 51 8 766 1 E-03, - 7. 9302 1 1 9 54E- 1 5 

i • 200000000E+00 , 0 . OOOOOOOOOE+OO : 3 . 886428 43 7E-03, - 5 . 570 740 690 E- 1 5 
8 • 103408 722E-03, -9 . 68 1 39 42 1 3E- 1 5 : - 4. 2 1 698028 6 E-0 3, 4. 1 1 06 53 523 E- 1 5 

1 .300000000E+00, 0 • OOOOOOOOOE + OO : 5 . 0 1 22 76 1 34E- 0 3, - 1 . 8 4 1 70 5356E- 1 4 
1 .041393 414E-0 2,-5.366821725E-14 : - 5 . 40 1 6 580 1 3E-0 3, 3 . 5251 1 63 70F- 1 4 

1 .400000000E + 00, 0 • OOOOOOOOOE + OO : 1 .02082 1 592E-0 1 , 4 . 9 4202 1 580E- 11 
-6.230056143E-03,-7.275944016E-12 : 1 . 083 1 22 1 53E- 0 1 , 5 . 669 6 1598 1 E- 1 1 

1 . 500Q00000E + 00, 0 • OOOOOOOOOE+OO : 8 . 8 508 5 458 6E + 0 2, 4.203559 133E-0 7 
1 .910474609E+02, 1 . 06 599802 7E-0 7 : 6 .9403799 77E+02, 3 . 1 3 756 1 1 0 6E- 0 7 


INPUT: IG0T0 
?4- 

Du4 9 


typical, run 

a. find region of root 




4 I NF U 1' : IT fp 


6 I i\ ! Ftj > I r, i N 1 £ i . El r; 1 r u. i\'El ;\ 1 bi . Chr i\! 1 L 
i 1 • A d j 0 • U 0.1* 0 E — 3 . 0 • 0 0 . 2 0 . 1 • 0 E - 6 


b. find root precisely 


1 Hxt’. F. ~ D i M Erl 5 I bi\'AL BAF F’L ED CHAM h EH 


k 1 = 1*00 Z 0 = 0 • 0 0 L 1 = 2 . 1 0 Z2= 2 • 8 0 

BETA0P= 0*00 0.00 BETA2L= 0-00 0.00 

l r\ N C:i £ iv F 1 A L M 0 1) £ i\' U M BEH 1 

M TO V i’j C = 5 MR 0 \/MC= 3 

wi"UVMC= o NRUVMC= 3 

J I T = 1 0 


El GiX ;< 1 
C H A K F N 

1 • 4b 0 0 0 0 0 0 0 £ + 0 0 . 0*0 0 0 0 0 0 0 0 0 E + 0 0 
2 . 4 6 6 63269 is £-0 1 . 3* 1 72 6 3 bO 1 4 E - 1 0 


1 . 4 4 3 7 0410 5 E+00. 

0 .00 0000000 E + 00 

• 

• 

5 .9220804942 + 00. 

2.934867061 

IF- 09 

6 .019 9 o 4 2 2 7 E ~ 0 2 . 

1 . 6 3200 39 72 E- 1 0 

• 

• 

5 . 8 568809 52 E + 00. 

2 . 769666464 

E - 0 9 

1 » A A 3 8 0 0 4 0 o E + 0 0 . 

0 .OOOOOOOOOE + OO 

J 

4 . 8 49 42620 5E + 00 . 

2.40 71 40 180 

E - 0 9 

1 » 0 7 3 4 7 9 9 8 7 E - 0 S . 

1 . 1 40 781 782 E- 10 

* 

4.83869 1 405E+00. 

2 . 293062002 

E - 0 9 

1 • 4 A 8 9 71 4 0 6 E + 0 0 . 

0 .000000000E + 00 

: 

4.622473622E+00. 

2 • 29 5354934 

[.; - o 9 

4*9 5886 1 1 0 0 E-04. 

1 .0 385181 76 E- 1 0 

; 

4*621977 736 E+00. 

2. 19 1 5031 1 7 

£-09 

1 * 4 4 2 9 4 4 8 6 6 E + 0 0 . 

0 .OOOOOOOOOE+OO 

2 

4.61 1225931 E+00. 

2.2898 1 3620 

£-09 

1 .2800323 i 2 E “06. 

1 .0 33 508 791 E- 10 

* 

4.61 122471 1 E+00. 

2. 186462741 

E - 0 9 

1 • 4429 4449 92 + 00. 

•0 .OOOOOOOOOE+OO 

: 

4.6111981 52E+00. 

2 • 289 799 934 

E-09 

2 *9103 8 3 0 4 6 £ - 1 1 . 

1 • 0 3349 6426 E- 10 

: 

4.6111981 52 E+00 . 

2. 18 6450291 

£-09 


I HP UT : I GOTO 


S UM 1 
SUMS 

: 8 . 7 A 8 b 3 7 b 6 0 S + 00. 4.32162 3 2 9 9 E - 0 9 
: 8.301 3 7 A 2 9 0 E + 0 0. 4.0 0 4 3 5 9 7 9 8 E - 0 9 
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4 INPUT: I TYPE 
?0 

5 I NPUT : El R ! ST j El R 1 SP > NEI R1 S j I DETAL 
?1 .4429 4 449 9., 0-00., 0-00 .>0.00* b 1 

INPUT: I PBSUM * IPV8UM* 1PAS0M, IPVAOM 
?Q * 0 * 0 * 0 

I NP UT : MP BSUM * MPV8UM* MPASOM * MP VAOM 
? 1 * 0 * bO 


THREE-DIMENSIONAL BAFFLED CHAMBER 


MUBAR= 3 

Rl= 1.00 Z0= 0*00 Zl= 2-10 Z2= 2-80 

BETAOP= 0-00 0.00 BETA2L= 0.00 0.00 

TANGENTIAL MODE NUMBER 1 . 

MTOVMC- 5 MROVMC= 3 

NTUVMC= 5 NRUVMC= 3 

JI T= 1 0 


EIGNR1 SUM 1 

CHARFN SUM 2 

1 .442944499E + 00* 0 .OOOOOOOOOE + OO : 4 • 6 1 1 1 9 7980 E + 00 * 2 . 289 7998 49 E- 09 
-7.68341 1241E-09, 1 . 0 33 49 63 48 E~ 1 0 : 4 . 6 1 1 1 9 7988E+00* 2 . 18645021 4E-09 


check pressure match 
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BSUM 



I NPUT : NRPNT 

s TVALj NTPNT* RVAL* TST0P ?20* 1 *0E-8> 40* 1 .0* 1 *05 

R- VAR I ATI 0N COMPARTMENT 1 TVAL= 0*0000*PI 

3 .208488334E-01 

X 

0*000000000E+00 


3 « 68066 788 1 E - 0 1 

5. 000000 000E-02 


4* 522753703E-01 

1 .OOOOOOOOOE-O 1 


5.514674265E-01 

1 • 500000000E-01 


6 . 5424671 19E-01 

2 .OOOOOOOOOE-O 1 


7 . 53840 1985E-01 

2. 500000000E-01 


8 • 461021 757E-01 

3 .OOOOOOOOOE-O 1 


9.281 785572E-01 

3. 500000000E-01 


9 .9802Q6354E-01 

4. OOOOOOOOOE-O 1 


1 .054651052 E+00 

4. 500000000E-0 1 

p( M )( £ )( r ,0.0,z 1 ) 

1 -098659082E+00 

5« OOOOOOOOOE-O 1 

1 . 1 323329 55E+00 

5. 500000000E-01 


1.15911 58 10E+00 

6 • OOOOOOOOOE-O 1 


1 * 1 82539420E+00 

6.500000000E-01 


1.2051 18913 E +00 

7 . OOOOOOOOOE-O 1 


1 -22771 6042E+00 

7. 500000000E-0 1 


1 « 249 60 7361 E + 00 

8 .OOOOOOOOOE-O 1 


1 .269 1 1 51 44E+00 

8. 506000000E-01 


1 .284408367E+00 

9 .0000000 00E-01 


1 . 294086334E+00 

9.500000000E-01 


1 .29738331 6E+00 

1 .OOOOOOOOOE + OO 


THETA- VARIATION 

RVAL= 1 .OOOOOOOOOE+OO 

1 *29738331 6E+00 

U 

0*0000*PI 

1 

1 « 282450983E + 00 

0*0 500*P I 

1 

1 .23022S410E+00 

o. 1000*PI 

1 

1 . 1 32 53826 5E+00 

0 * 1 500*PI 

1 

1 .006 7101 8 7E+00 

0*2000*PI 

1 

9 .0 1 399 5673E-01 

0 * 2500*P I 

1 

8 .654279342E-01 

0 • 3000+P I 

1 

8.991 267796E-01 

0 • 3500*P I 

1 

9 .31 3336833E-01 

0 *4000*PI 

1 

8 • 524209203E-01 

0 . 4500*PI 

1 f n\f 0\ J 

5.90751 9324E-01 

0 • 5000*P I 


1 .802798582E-01 

0 * 5500*PI 

1 

-2. 351 8 90570 E- 01 

0 • 6000*P I 

1 

-4.741978769E-01 

0 *6500*PI 

1 

-8 .6734751 64E-0 1 

0. 7000*PI 

2 

- 1.1 44832867E+00 

0 • 7500*PI 

2 

- 1 . 49633349 4E + 00 

0 *8000*PI 

2 

- 1 • 747986765E+00 

0 *8500*PI 

2 

-1 .827863061 E+00 

0 • 9000*P I 

2 

- 1 • 79 5720557E + 00 

0 .9500*PI 

2 

- 1 *7663573 59 E+00 

1 *0000*PI 

2 

- 1 . 7957205 5 7E+00 

1 *0500*PI 

2 

-1 *827863061 E+00 

1 * 1000*PI 

2 

I NPUT : NRPNT > TVALj NTPNT, RVAL j TST0P ?0* 0 *0* 0> 0 * 0*0*0 
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I i\IP U 1 I NRi 5 NT > T V AL j NTP NT* H’V AL t T S Fur ? 20 * 1 • 0 h- 8 > > i • ■ 1 j 1 • 0 '> 


~ V A K I A l I b M i v l A I >\ ; C H A ft B E K 


Tv Ai 


1 • 3 72 3 3 60 2 2 E- 12 
2*1911 2 ft 7 3 3 E - 0 1 
ft • 0 6 1 5 6 6 6 3 7 E - 0 1 
5*60g279405E-0l 
6 • 8 529 78 ft 6 7 £ - 0 1 
7 • ft 7 79 6 7 0 2 3 E - 0 1 
6 * 7 40 8 10 18 1 E-CJ1 
9 • ft 6 6 2 9 8 4 0 ft E - 0 1 
1 • 0 1 345 39 £! E + 00 
1 • U 6 6 7 1 59 2 1 E + 0 0 
1 • 1 1 ft 0 5 0 8 09 E+00 
1 » 1 49 9 3 30 1 ft E + 0 0 
1 • 1 78 31 860 1 £+00 
1 • 20321 4786 8 + 00 
1 *227 b 6 6 7 6 0 E + 0 0 
1 • 2 3 3 1 6 3 0 9 2 E + 0 0 
1 • 28 j 49 3 16 3E + 0Q 
1 • 3 1 6 58 90 78 E+00 
1 • 3 ft ft 1 1 50 76 E + 00 
1 • 3 6 3 1 6 9 3 2 ft E + 0 0 
1 * 3 699 771 8 5 E + 00 


r 

0 *000000000 E + 00 
3 * 0 0 0 0 0 0 0 0 0 E - 0 2 
1 *0000 0000 OE-O 1 

1 * 5 0 0 0 00000 E - 0 1 
2* 000000000 E- 01 
2 • 30000 0 0 0 0 E - 0 1 
3* 0000 0000 0E- 01 
3* 50000000 OE-O 1 
ft* 00 00 00 00 OE-O 1 
ft *30 0 0 i ) 0 0 0 0 E - 0 1 
3* 000000000 £-0 1 
3* 30 00 0000 OE-O 1 

6 *0000000 0 0 E - 0 1 
6 • 30 0 0 0 0 0 0 OE-O 1 

7 * 0 0 0 0 0 00 0 0 E - 0 1 

7 • 50 0 0 0 0 0 0 0 E - 0 1 

8 * 00000000 OE-u 1 

8 * 500000000E-0 1 
9*00 0 0 0 0 0 0 0 E - 0 1 

9 * 300000000 E-0 1 
1 .OOOOOOOOOE+OO 


0 • 00 00 sfCp I 


p(0)( e )(r JI O.O,z 1 ) 


THETA- 9ARI AT I ON 

1 *3699771 8 5E+00 
1 * 3 434371 19 E + 00 
1 * 2 72900 338 E + 00 
1 * 181 /Q 999 6 E + 00 
1 *09 729 620 3 E + 00 
1 • 0 38 78 58 43E + 00 
1 *007379 760 E+00 
9 *o3ft8oft356E-01 
9 * 3 2 0 6 3 8 4 7 5 £ - 0 1 
8 * 1 4443 5 752E-0 1 
6 *02441 5387E-01 
2*895868696E-Q1 
~ 1 *04 /23o489E~01 
-5.397050468E-01 
-9*6471 1 4272E-01 

- 1 *3336436 18 E + 00 

- 1 *61 70544 40 E+00 
-1 *80772551 2E+00 
-1 * 9 1 SO 1021 OE + OO 

- 1 -970377663E + 00 
-1 -985144713E+00 

- 1 *9 70377663E + 00 
-1 *918010 210 E + 00 


RVAL = 

e 

0 *00 00 ftp I 
0 *050 Oft PI 
0 * 1000 ftp I 
0 * 1 bOOftPI 
0 • 2000 ftp I 
0 *2500 ftP I 
0 *3000 ftp I 
0 *3 500ftPI 
0 * 40 00 ftp I 
0 * 4 500 ftP I 
0 * SOOOftPI 
0* 5500ftPI 
0 * 60 00 ftp I 
0*6 5Q0ftp I 
0 • 7000ftp I 
0 * 7 500 ftP I 
0* SOOOftPI 
0 *8500ftPI 
0 • 9000 ftp I 

0 *9500*PI 
1 .OOOOftPI 
1 *0500ftPI 

1 * lOOOftPI 


•OOOOOOOOOE+OO 


p(O)(C)( 1#Oj0j2i ^ 


I NPUT : NRPNTj TVAL., NTPiMT* RVAL* TST0P ?0*0.0,0*0*0*0.6 
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NOMENCLATURE FOR APPENDICES C AND D 


A s series expansion coefficient - pressure gradient 

a s series expansion coefficient - pressure 

B s series expansion coefficient - pressure gradient 

b s series expansion coefficient - pressure 

c s sound speed 

F e axial Green’s function term 

G = Green’s function 

H a axial influence coefficient 

I = integral coefficient 

J s Bessel function 

k = eigenvalue 

m = summation index 

n s summation index 

p s pressure 

r = radial coordinate - 3D chamber 
x =. axial coordinate - 2D chamber 
y s transverse coordinate - 2D chamber 
z s axial coordinate - 3D chamber 

a = rigid wall eigenvalue 
P = specific acoustic admittance 

r = Gamma function 
6 = phase angle 

6 = angular coordinate - 3D chamber 
A = normalization parameter 

UJ s eigenfunction 

v-5k 



Subscripts 


0 = nozzle end 

1 = interface 

2 s injector end 
L = left 

m 5 summation index 
n = summation index 
R = right 

r = radial component (3D) 
x = axial component (3D) 
y = transverse component (3®) 
z = axial component (3D) 

p 5 compartment number 
p = number of compartments 

Superscripts 

(0) = main chamber 

(p) s p th compartment 
( C ) = iteration number 



EQUIVALENCE OF SYMBOLS USED IN MAIN 


TEST AND APPENDICES C AND D 


MAIN TEXT 

APPENDICES 

c 

c 

G 

G 

k 

k 

k 

. q 

(p) 

a 

X,n j 

k 

m 

* (0) 

x,m y 

2 

X 2 ~ X 1 > Z 2 ~ Z 

L 

x -x . z, -z 

1 O 7 1 1 

zero 

x . z 
o 7 0 

m 

m -1 

y 

hitt/W 

1(0) 

2 y>"v 

p a 

p(0) 

P bM- 

p w 

q 

n -1 

y 

qn/v 

(p) 

a 

y,Hy 
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MAIN TEXT 


R 


w 


W 


zero 


N 


APPENDICES 


Vl“^ 


y_-y 
n o 


2 - 


o+ 


D-5T 



DISTRIBUTION LIST 


Dr. R. J. Priem MS 500-209 
NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland., Ohio 44135 (2) 

Norman T. Musial 
NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland , Ohio 44135 

Library 

NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland^ Ohio 44135 (2) 

Report Control Office 
NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland; Ohio 44135 

Brooklyn Polytechnic Institute 
Attn: V. D. Agosta 

Long Island Graduate Center 
Route 110 

Farmingdale; New York 11735 

Chemical Propulsion Information Agency 

Johns Hopkins Uni vers ity/APL 

Attn: T. W. Christian 

8621 Georgia Avenue 

Silver Spring; Maryland 20910 

NASA 

Lewis Research Center 
Attn: E. W. Conrad; MS 500-204 

21000 Brookpark Road 
Cleveland; Ohio 44135 

North American Rockwell Corporation 

Rocketdyne Division 

Attn: T. A. Coultas; D /578 BA28 

6633 Canoga Avenue 

Canoga Park; California 91304 

National Technical Information Service 
Springfield; Virginia 22151 
(40 Copies) 

NASA Representative 
NASA Scientific and Technical 
Information Facility 
P.0. Box 33 

College Park; Maryland 20740 


Aerospace Corporation 
Attn: 0. W. Dykema 

Post Office Box 9508.5 
Los Angeles; California 90045 

Ohio State University 
Department of Aeronautical and 
Astronautical Engineering 
Attn: R. Edse 

Columbus; Ohio 43210 

TRW Systems 

Attn: C-. W. Elverum 

One Space Park 

Redondo Beach; California 90278 

Bell Aerospace Company 
Attn: T. F. Ferger 

Post Office Box 1 
Mail Zone J- 8 l 
Buffalo; New York 14205 

Pratt & Whitney Airci raft 

Florida Research & Development Center 

Attn: G. D. Garrison 

Post Office Box 710 

West Palm Beach; Florida 33402 

NASA 

Lewis Research Center 
Attn: L. Gordon; MS 500-209 

21000 Brookpark Road 
Cleveland; Ohio 44135 

Purdue University 

School of Mechanical Engineering 

Attn : R. Goulard 

Lafayette; Indiana 47907 

Air Force Office of Scientific Research 

Chief Propulsion Division 

Attn: Lt. Col. R. W. Haf frier (NAE) 

1400 Wilson Boulevard 
Arlington, Virginia 22209 

Pennsylvania State University 
Mechanical Engineering Department 
Attn: G. M. Faeth 

207 Mechanical Engineering Bldg. 
University Park, Pennsylvania 16802 


(2 Copies with Document Release Authorization Form) 



; i 


University of Illinois 
Aeronautics/ Astronautic Engineering 

Department 

Attn: R. A. Strehlow 

Transportation Building, Room 101 

Urbana, Illinois 6l801 

NASA 

Manned Spacecraft Center 
Attn: J. G. Thibadaux 

Houston, Texas 77058 

Massachusetts Institute of Technology 
Department of Mechanical Engineering 

Attn: T. Y. Toong 

17 Massachusetts Avenue 
Cambridge, Massachusetts 02139 

Illinois Institute of Technology 

Attn: T. P. Torda 

Room 200 M. H. 

3300 S' Federal Street 
Chicago, Illinois 60616 

AFKPL 

Attn: R. R, Weiss 

Edwards, California 93523 

LAS' Amy Missile Command 
AMSMI-RKL, Attn: W. W. Wharton 

Redstone Arsenal, Alabama 35808 

University of California 
Aerospace Engineering Department 

Attn: F. A. Williams 

Post Office Box 109 
La Jo 11a, California 92037 

Georgia Institute of Technology 

Aerospace School 
Attn: B. T. Zinn 

Atlanta, Georgia 30332 

Marshall Industries 
Dynamic Science Division 
Attn: L. Zung 

2400 MLchelson Drive 
Irvine, California 92664 

Mr. Donald K. Dahlene 
U. S. Army Missile Command 
Research, Development , Engineering 
and Missile Systems Laboratory 

Attn ; AMSMI-RK 

Redstone Arsenal, Alabama 35809 


TISIA 

Defense Documentation Center 
Cameron Station 
Building 5 
5010 Duke Street 
Alexandria, Virginia 22314 

Office of Assistant Director 
(Chemical Technician) 

Office of the Director of Defense 
Research and Engineering 
Washington, D. C. 20301 

D. E. Mock 

Advanced Research Projects Agency 
Washington, D. C. 20525 

Dr. H. K. Doetsch 

Arnold Engineering Development Center 
Air Force Systems Command 
Tullahoma, Tennessee 37389 

Library 

Air Force Rocket Propulsion Laboratory 
(RPR) 

Edwards, California 
Library 

Bureau of Naval Weapons 
Department of the Navy 
Washington, D. C. 

Library 

Director (Code 6l8o) 

U. S. Naval Research Laboratory 
Washington, D. C. 20390 

APRP (Library) 

Air Force Aero Propulsion Laboratory 
Research and Technology Division 
Air Force Systems Command 
United States Air Force 
Wright-Patterson AFB, Ohio 45433 

Technical Information Department 
Aeronutronic Division of Philco Ford 
Corporation 
Ford Road 

Newport Beach, California 92663 

Library- Do cuments 
Aerospace Corporation 
2400 E. El Segundo Boulevard 
Los Angeles, California 90045 


Princeton University 
James Forrestal Campus Library 
Attn: D. Harrje 

Post Office Box 710 
Princeton, New Jersey 08540 

U. S. Naval Weapons Center 
Attn: T. Inouye, Code 4581 

China Lake, California 93555 

Office of Naval Research 
Navy Department 
Attn: R. D. Jackel, 473 

Washington, D. C. 20360 

Air Force Aero Propulsion Laboratory 
Attn: APTC Lt. M. Johnson 

Wright Patterson AFB, Ohio 45433 

Naval Underwater Systems Center 
Energy Conversion Department 
Attn: Dr. R. S. Lazar, Code TB 131 

Newport, Rhode Island 02840 

NASA 

Langley Research Center 
Attn: R. S. Levine, MS 213 

Hampton, Virginia 23365 

Aerojet General Corporation 
Attn: J. M. McBride 

Post Office Box 15847 
Sacramento, California 95809 

Colorado State University 
Mechanical Engineering Department 
Attn: C. E. Mitchell 

Fort Collins, Colorado 80521 

University of Wisconsin 
Mechanical Engineering Department 
Attn: P. S. iYfy-ers 

1513 University Avenue 
Madison, Wisconsin 53708 

North American Rockwell Corporation 
Rocketdyne Division 
Attn: J. A. Nestlerode 

Ac46 D/596-124 

6633 Canoga Avenue 

Canoga Park, California 91304 


University of Michigan 
Aerospace Engineering 
Attn: J. A. Nicholls 

Ann Arbor, Michigan 4810 4 

Tulane University 
Attn: J. C. O'Hara 

6823 St. Charles Avenue 
New Orleans, Louisiana 70118 

University of California 
Department of Chemical Engineering 
Attn: A. K. Oppenheim 

6l6l Etcheverry Hall 
Berkeley, California 94720 

Army Ballistics Laboratories 
Attn: J. R. Osborn 

Aberdeen Proving Ground, Maryland 21005 

Sacramento State College 
School of Engineering 
Attn: F. H. Reardon 

6000 J. Street 

Sacramento, California 958X9 

Purdue University 

School of Mechanical Engineering 

Attn: B. A. Reese 

Lafayette, Indiana 47907 

NASA 

George C. Marshall Space Flight Center 
Attn: R. J. Richmond, SNE-ASTN-PP 

Huntsville, Alabama 35812 

Jet Propulsion Laboratory 
California Institute of Technology 
Attn : J. H. Rupe 

4800 Oak Grove Drive 
Pasadena, California 91103 

University of California 
Mechanical Engineering Thermal Systems 
Attn: Prof. R. Sawyer 

Berkeley, California 94720 

ARL (ARC) 

Attn : K. Scheller 

Wright Patterson AFB, Ohio 45433 



Library 

Bell Aerosy stems, Inc. 

Box I 

Buffalo, New York 14205 

Report Library, Room 6A 
Battelle Memorial Institute 

505 King Avenue 
Columbus , Ohio 43201 

D. Suichu 

General Electric Company 

Flight Propulsion Laboratory Department 

Cincinnati, Ohio 45215 

Library 

Ling- Temco - Vo ugh t Corporation 

Post Office Box 5907 
Dallas, Texas 75222 

Marquardt Corporation 
16555 Saticoy Street 
Box £013 - South Annex 
Van Nuys, California 91409 

P. F. Winternitz 
New York University 
University Heights 
New York, New York 

I. Forsten 
Picatlnny Arsenal 
Dover, New Jersey 078 OI 

R. Stiff 

Propulsion Division 
Aerojet-General Corporation 
Post Office Box 15847 

Sacramento , California 95803 

Library, Department 598 - 306 
Rocketdyne Division of Rockwell 
North American Rockwell Inc. 

6633 Caaoga Avenue 

Canoga Park, California 91304 

Library 

Stanford Research Institute 

333 Ravenswood Avenue 
Menlo Park, California 94025 


Library 

Susquehanna Corporation 
Atlantic Research Division 
Shirley Highway and Edsall Road 
Alexandria, Virginia 22314 

STL Tech. Lib. Doc. Acquisitions 
TRW System Group 
1 Space Park 

Redondo Beach, California 90278 

Dr. David Altman 
United Aircraft Corporation 
United Technology Center 
Post Office Box 358 
Sunnyvale, California 94088 

Library 

United Aircraft Corporation 

Pratt and Whitney Division 

Florida Research and Development Center 

Post Office Box 2691 

West Palm Beach, Florida 33^02 

Library 

Air Force Rocket Propulsion 
Laboratory (RPM) 

Edwards, California 93523 


Allan Hribar, Assistant Professor 
Post Office Box 5014 
Tennessee Technological University 
Cookeville, Tennessee 38501 

NASA 

Lewis Research Center 
Attn: E. 0. Bourke MS 500-209 

21000 Brookpark Road 
Cleveland, Ohio 44135 

NASA 

Lewis Research Center 
Attn: D. L. Nored 500-203 

21000 Brookpark Road 
Cleveland, Ohio 44135 

NASA 

Lewis Research Center 
Attn: Contracting Officer, 500-313 

21000 Brookpark Road 
Cleveland, Ohio 44135 



